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SHORT-TERM ASYMPTOTICS FOR THE IMPLIED VOLATILITY SKEW UNDER A 
STOCHASTIC VOLATILITY MODEL WITH LEVY JUMPS 

JOSE E. FIGUEROA-LOPEZ AND SVEINN OLAFSSON 


Abstract. The implied volatility skew has received relatively little attention in the literature on short-term asymptotics 
for financial models with jumps, despite its importance in model selection and calibration. We rectify this by providing 
high-order asymptotic expansions for the at-the-money implied volatility skew, under a rich class of stochastic volatility 
models with independent stable-like jumps of infinite variation. The case of a pure-jump stable-like Levy model is also 
considered under the minimal possible conditions for the resulting expansion to be well defined. Unlike recent results for 
“near-the-money” option prices and implied volatility, the results herein aid in understanding how the implied volatility 
smile near expiry is affected by important features of the continuous component, such as the leverage and vol-of-vol 
parameters. As intermediary results we obtain high-order expansions for at-the-money digital call option prices, which 
furthermore allow us to infer analogous results for the delta of at-the-money options. Simulation results indicate that 
our asymptotic expansions give good fits for options with maturities up to one month, underpinning their relevance 
in practical applications, and an analysis of the implied volatility skew in recent S&P500 options data shows it to be 
consistent with the infinite variation jump component of our models. 
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1. Introduction 

1.1. Motivation. Since the emergence of the Black-Scholes option pricing model, there have been two main themes 
in the evolution of hnancial models: jumps and stochastic volatility. The features of both model components have 
their merits, and combining jumps in returns and stochastic volatility is a viable way to calibrate the implied volatility 
surface across strikes and maturities. However, the increased generality comes at a cost as models that combine jumps 
and stochastic volatility are often highly complicated to calibrate and implement. An important line of research has 
therefore been to examine various extreme regions of the volatility surface and attempt to understand how different 
model features and their associated parameters affect the behavior of option prices (see [4], [5], [15], [17], [18], [2T], 
[22] . [31], [39], and references therein). 

In particular, a number of recent papers have shed light on the short-term asymptotic behavior of option prices and 
implied volatility, and revealed that those quantities exhibit markedly different behavior from one model setting to 
the next. For instance, while there typically exists a limiting implied volatility smile in continuous models, the smile 
exhibits explosive behavior in models with jumps. This phenomenon is due to the much slower decay of out-of-the- 
money (OTM) option prices in the presence of jumps, which causes the implied volatility to blow up off-the-money, 
while converging to the spot volatility at-the-money (ATM) (see, e.g., [T7], [5T], and [SH])- This feature of jump models 
is actually desirable as it enables them to reproduce the pronounced smiles and skews observed at short maturities. 

Another stylized empirical fact is that as time-to-maturity decreases, the liquid strike prices become increasingly 
concentrated around the ATM strike (see, e.g., [3T], and Figure herein). ATM options are therefore of particular 
importance when it comes to short-term asymptotics. The leading order term of ATM option prices has been derived 
for various models containing jumps (|32j. |39| 1. but, unfortunately, those approximations are known to require unreal¬ 
istically small maturities to attain satisfactory accuracy, and provide limited information about key model parameters. 
These drawbacks have motivated the search for higher order asymptotics and, in the recent paper m, such expansions 
are derived for a general class of stochastic volatility models with Levy jumps, in an asymptotic regime where time- 
to-maturity and log-moneyness become small, which, as explained above, is of particular practical importance. Still, 
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certain shortcomings prevail, in particular when it comes to understanding the effect of the stochastic volatility com¬ 
ponent on the short-term volatility smile, and its interaction with the jump component. For instance, the asymptotic 
expansions developed in [ 21 ] depend on the volatility process only through the spot volatility. 

In the present work, we alleviate the issues mentioned in the previous paragraph by analyzing the short-term behavior 
of the ATM implied volatility skevQ which turns out to depend on key parameters of the underlying volatility process 
such as the leverage and vol-of-vol parameters, and gives a more accurate picture of the short-term behavior of the 
implied volatility smile. The ATM skew has received comparatively little attention in the literature, despite the fact 
that it is actively monitored in practice by traders and analysts (cf. [30] 1 . stemming in part from its rich informational 
content, and considerable empirical support has in fact been provided for its significance in predicting future equity 
returns, and as an indicator of the risk of large negative jumps (see, e.g., |40j . [41] . for individual stock options, and 
| 6 ], [34], for index options). 

The ATM skew is also highly relevant in model selection and when calibrating models to observed option prices (cf. 
[23l Ch. 5]), especially in FX markets, where it is standard to effectively quote directly on smile slope and convexity 
(cf. [3 Sec. 2]). Moreover, it is also widely believed that the implied volatility skew and convexity reflect the skewness 
and kurtosis in the underlying risk-neutral distribution and, thus, a number of researchers have attempted to relate 
the smile skew and convexity to moments of the risk-neutral distribution (see, e.g., m. Ba, m, m)- 


Finally, another important reason to study the skew is its connection with the delta of options, which is of paramount 
importance in the trading and hedging of options. More specifically, the same key quantities needed to derive the 


asymptotic behavior of the skew, can be used to derive short-term asymptotics for the delta of options (see Section 1.3 
below). 


1.2. Literature Review. The literature on the short-term implied volatility skew in the presence of jumps is somewhat 
limited, but in recent years the leading order term has been obtained, albeit under somewhat restrictive assumptions. 
For models of jump-diffusive nature, several sources show that the skew converges to a nonzero value as time-to-maturity 
tends to zero (see, e.g., [ 1 ], [TS], [II])- In such models, jumps are infrequent and can be interpreted as the occurrences 
of rare events, but substantial empirical evidence supports models with infinite jump activity. In particular, S&P500 
index options are used in [53 to reach the conclusion that adding Poisson jumps to a stochastic volatility model is not 
sufficient to account for the implied volatility skew at short maturities. 


For models with infinite jump activity, a recent result m shows that the skew is of order t for bounded variation 
Levy processes, as well as a few specific infinite variation cases with Blumenthal-Getoor inde?!^ equal to one, such as 
the Normal Inverse Gaussian (NIG) and Meixner processes. This growth rate is in fact the fastest possible one in the 
absence of arbitrage (cf. [28[ Sec. 3]), but in [15] the authors also show that for certain Levy models with a Brownian 


component, the skew can explode at a rate slower than —1/2. This can be viewed as a special case of Corollary 4.3 in 
Section [^herein, and the same can be said about Proposition 8.5 in a recent survey paper [3, which can be interpreted 
as a skew approximation for tempered stable Levy processes as defined in |12j . but under some extremely restrictive 
assumptions on the model parameters. 


A somewhat different approach is adopted in m, which belongs to the stream of literature attempting to relate 
features of the implied volatility smile to properties of the risk-neutral distribution of the underlying. Concretely, the 
authors derive, partly by heuristic arguments, the following expansion for near-the-money implied volatility 
parameterized in terms of the scaled moneyness k = {K — Sq)/{SQ( 7^/i), where K and t are the option’s strike and 
time-to-maturity, Sq is the spot price of the underlying, and cr is a measure of the “overall volatility level”: 


a{k,t) = a [at + Ptk + "ftk"^ + 0{k^)). 


(l.I) 


Explicit expressions for the coefficients at, I3t, and 74 are also proposed (see Eq. (2) therein). In particular, the form 
of the skew coefficient fit suggests, at least heuristically, the following expression for the ATM skew. 


da (k, t) 


dk 


= a 




2-P(^i>0) 


+ o{t 2 ) ^ t — 0 , 


( 1 . 2 ) 


which agrees with formula (1.3) below, but can at best be used to obtain the leading order term of the skew. More 
importantly, a mathematically sound justification is needed when passing from the asymptotic expansion ( 1 . 1 ) for 


^Practitioners commonly use the terms “skew” and “implied volatility skew” for the ATM slope of the implied volatility curve for a given 
expiration date (see, e.g., m)- We will use the terms interchangeably. 

^For a Levy process X with Levy measure u, the Blumenthal-Getoor index is defined as inffp > 0 : \x\Pu(dx) < oo}. 
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d{R,t) to an analogous asymptotic expansion for its derivative of the form (1.2). In a certain sense, the results herein 
therefore formalize and extend the heuristic approach for the skew in |13j . Another important contribution of the 
approach in [13] is that it attempts to explain the general shape of volatility smiles and, in particular, how the smile 
skew, j3tj is related to the skewness of the distribution of the underlying. To this end, the authors use S&P500 index 
time series to argue that there cannot be a simple relation between the skew of the smile and the skewness. 


1.3. Overview of New Results. It is important to stress that most of the aforementioned results still suffer from the 
same shortcomings as the corresponding results for ATM option prices, in that their domain of validity is extremely 
small. This, in part, stems from the focus being on obtaining the leading order term, which can be done in some 
generality, but in return tends to only depend on the most general model parameters. A significant contribution of 
the present work is therefore to provide accurate higher order expansions for the implied volatility skew, under a class 
of models that goes beyond the homogeneous Levy framework by combining stochastic volatility and jumps with high 
activity. Empirical evidence generally supports the need for such models, not only to calibrate the implied volatility 
surface, but also to generate realistic future dynamics of implied volatility, in order to give reasonable prices for exotic 
derivatives (see, e.g., 0)- 


Throughout, we assume that the risk-free interest rate r and the dividend yield <5 of the underlying are 0, and that 
the price process S := {St)t>o of the underlying asset is a P-martingale. We denote the implied volatility of an option 
by a(K,t), where k := \og{K/So) is the log-moneyness and t is the time-to-maturity. For simplicity, the ATM implied 
volatility (t( 0,<) is denoted by a{t). Let us start by recalling some basic relationships that are fundamental to our 
approach (see Section]^ for details). First, under some mild conditions the ATM skew satisfies, 

> 2 


da {K,t) 


dn 


271/1 


K=0 




Hthfi 

2 v^ 


O 


((o-(t)Vt)' 


1 -f 


'r{t)Viy 


o 


((CT(t)v^)^) 


(1.3) 


which enables us to separate the problem of studying the asymptotic behavior of the skew into finding the asymptotics 
of two quantities that are important in their own right: (i) the ATM implied volatility a{t)^ and (ii) the ATM digital 
call option price F {St > Sq). Interestingly enough, there is also a close connection between P(S't > S'o) and the delta 
of ATM options, i.e. the sensitivity of the ATM option price C{So,t), with respect to the spot price of the underlying, 
Sq. Concretely, we show that 


A{t) := ^£^2^ = ^c{So,t)+F{St > So), 

ODq Oq 


(1.4) 


so the asymptotic results for the ATM option price C'(iS'o, t) (cf. [H]), together with the present results for the transition 
probability P(5't > So), can be used to obtain short-term asymptotic expansions for the delta of ATM call options. 

Both the ATM implied volatility, a{t), and the corresponding option price, C{So,t), have received considerable 
attention in the literature, and it is well documented that their short-term behavior is strongly tied to various pathwise 
properties of the log-returns process. For example, including a continuous component can significantly change the 
properties of pure-jump models, and the type of jump component can also have a drastic effect. The same is true 
when it comes to the ATM volatility skew, and we proceed to explain the different novelties of our work by separately 
analyzing the two cases of interest (pure-jump and mixed), and then elaborating on the accuracy and applicability of 
our results in model selection and calibration. We also briefly consider the OTM volatility skew, which, much like the 
OTM volatility (cf. [TT]), can be analyzed in much more generality than the ATM skew (see Remark 4.5). 


1.3.1. Pure-jump exponential Levy model. In this paper, we consider tempered stable-like Levy processes, as introduced 
in m and [21] ■ with and without an independent continuous component. In this section we briefly describe our results 
related to the latter case. More concretely, we consider the model 

St:=Soe^\ (1.5) 

where X stands for a pure-jump Levy process with a Levy measure of the form 

^{dx) = c(-^yx\~^~^q{x)dx, ( 1 . 6 ) 

for some constants (7(1), (7(—1) S [0, oo) such that (7(1) -I- (7(—1) > 0, F S (0, 2), and a bounded function q : K\{0} — 
[0, oo) such that q{x) —>■ 1, as x —>■ 0. This framework includes most of the infinite activity Levy models used in 
practice, and its short-term behavior depends strongly on the index of jump activity Y, which coincides with the 
Blumenthal-Getoor (BG) index of the process. In what follows we impose the condition Y G (1, 2), which implies that 
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X has infinite variation. This restriction is supported by recent econometric studies of high-frequency financial data 
(see Remark 2.2 in m), and we will also argue in Section]^ using recent S&P500 options data, that it is the most 
relevant case since it gives the flexibility needed to calibrate the short-term implied volatility skew observed in practice. 


For models of the form (1.5), a second order short-term expansion for the ATM implied volatility, is given in 
Theorem 3.1 of [21j . under a minimal integrability condition on q around the origin. The key to studying the implied 
volatility skew is therefore the transition probability appearing in ( |1.3[ ), for which we have ¥ {St > Sq) = P(A^t > 0), 
but while a lot is known about P {Xt > xg) for nonzero xg (cf. [20]1. much less has been said about P (Xt > 0) for 
processes with infinite jump activity. The leading order term for bounded variation Levy models, as well as certain 
models with BG-index one (e.g., NIG, Meixner) is obtained in [TS], while, for tempered stable-like processes, 


P(Xt>0) — >P(Zi>0), 0, 


(1.7) 


where Zi is a strictly M-stable random variable under P. This limit is a consequence of the fact that Xt converges 
in distribution to Zi, as t —>■ 0 (cf. [36]), and cannot be extended to higher order terms. However, procedures similar 
to the ones used to derive near-the-money option price expansions in [18] and [2T] will allow us to get a closer look at 
the convergence. Concretely, the following novel higher order asymptotic expansion is obtained, 

n 

P{Xt>0)-P{Zi>0)='^dkt''(^-^')+et^+ft + o{t), t^O, ( 1 . 8 ) 

k=l 


where n := maxjfc > 3 : k{l — 1/M) < 1}. It is important to point out that the leading order term is for all 

Y G (1, 2). This can be compared to the expansion for ATM option prices given in Theorem 3.1 of [21], where the first 
and second order terms are of order and t, i.e. the convergence here of digital option prices is slower. 


With the expansion (1.8) at our disposal, we can use (1.3) to deduce an expansion for the ATM implied volatility 
skew, which turns out to exhibit explosive behavior in short time, but unlike jump-diffusion models where the skew 
is always bounded, and finite variation models where the skew is of order (cf. [13]), the order of the skew here 
actually depends on the index of jump activity, M, and ranges between 0 and —1/2. This provides an important 
model selection and calibration tool, and in Section [^ we will show that our result is in line with the short-term skew 
in S&P500 index options, while the same cannot be said about models where the dynamics are driven by a purely 
continuous model, nor a model with a finite variation jump component. Furthermore, important qualitative properties 
such as the sign of the skew can easily be recovered from the model parameters, and used to create simple parameter 
restrictions that can be used for calibration and model selection purposes (see Remark 3.4 for further details). 


1 . 3 . 2 . Exponential Levy model with stochastic volatility. Empirical work has generally supported the need for both 
jumps to reflect shorter maturity option prices, and stochastic volatility to calibrate the longer maturities where the 
smile effect of jump processes is limited (see, e.g., [531 Gh.5] and [5^). In order to incorporate a continuous component 
of diffusive type into the price dynamics, we consider the model 

St := Sge^*+^\ ( 1 . 9 ) 

where {Vt)t>g is a stochastic volatility process of the form 

dVt = yi{Yt)dt + aiYt){pdWt^ + y/l^dW^) , Mq = 0, 

dYt = a{Yt)dt + 'y{Yt)dWt, Id = yg, 

and {Wl)t>g and {Wt)t>o are independent standard Brownian motions, independent of the pure-jump Levy process 
X. This framework includes the most commonly used stochastic volatility models, such as the mean-reverting Heston 
and Stein-Stein models, and we remark that in such models it is generally believed that the leverage parameter p is 
responsible for generating asymmetric volatility smiles. A question of interest is therefore to what extent, and in what 
way, the leverage parameter contributes to the short-term skew in models with jumps. As in the pure-jump case, the 
key to studying the ATM skew is the transition probability P(S't > Sg)^ and for models of the form St = 
we have P(Xt -I- Id > 0) —>■ 1 / 2 , as t > 0. The literature is quite sparse beyond that, but in Section]^ we derive the 
following higher order expansion, 

P(A:t-lEt > 0) = --fiefs-fio(t^), t^O, (1.11) 
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where n := max{fc > 3 : k {1 — Yl2) < (3 — Y)/2}. Comparing this result to the short-term expansion for ATM call 
option prices given in Theorem 4.1 of [H], where the first and second order terms are of order and reveals 

that the convergence of ATM digital option prices is slower, as it was in the pure-jump case, unless C'(l) = C'(—1), in 
which case the summation term vanishes. 

Piecing together the results above now gives the following expansion for the ATM implied volatility skew. 


da {n,t) 


dn 


= -V^y^dkt(^ 2 )^ 2 - fy/2Trf + l^ai) ^ + o{t^ 0 , 

«=o ^ <^iyo) ^ 2 / 


( 1 . 12 ) 


where c := 7 — {pa'{yo)j{yo))/2, and 7 is a constant that depends only on the parameters of X. As in the pure-jump 
case, the order of the skew ranges between 0 and —1/2 depending on the index of jump activity, Y, but it is also 
observed that for a fixed value of T, the skew is less explosive than in the pure-jump case. In particular, in the 
symmetric case (7(1) = C'(—1), the skew actually converges to a nonzero value, —c/cr(i/o)j as in jump-diffusion models. 
This is to be expected since including a continuous component has a limited effect on OTM volatility, while raising the 
limiting ATM volatility from zero to the spot volatility a{yQ), effectively flattening the smile. The expression (1.12) 
also offers a significant improvement over existing results in that it depends on both the correlation coefficient p and 
the volatility of volatility, a'{yo)j{yo). In particular, this is in sharp contrast to the expansions for option prices and 
implied volatility in |2I| , where the impact of replacing the Brownian component by a stochastic volatility process was 
merely to replace the volatility of the Brownian component, a, by the spot volatility, a{yo)- 


Lastly, it is noteworthy that the skew-effects of the jump-component and the continuous component (i.e. the corre¬ 
lation coefficient p) turn out to be additive to the leading order, and the contribution to the implied volatility skew of 
a nonzero correlation p between the asset price and volatility can be quantified as 


1 p<^'iyo)i{yo) 

2 cr(?/o) 


( 1 + 0 ( 1 )), 


t ^ 0. 


Interestingly enough, this shows that to the leading order, the skew-effect of stochastic volatility is the same as in 
jump-diffusion models (see Ex. 7.1 in [3] for a comparable result for such finite activity models), but the effect of the 
jump-component is drastically different, and so is the interaction between the two model components. 


1.3.3. Accuracy of the asymptotic formulas/Empirical analysis. As explained above, the asymptotic expansions ob¬ 
tained herein can be used to infer several important features of the implied volatility smile. However, we also men¬ 
tioned that a common drawback of short-term approximations is that their domain of validity can be small, which 
could potentially limit their usefulness for practical work. To better assess this point, in Sectionj^we test the accuracy 
of our expansions using Monte Carlo simulation. Our results indicate that for the important class of tempered stable 
processes, they give good approximations for options with maturities up to one month, underpinning their practical 
relevance in essentially every major options market. In the second part of Sectionwe look at the short-term implied 
volatility skew in recent S&P500 option prices, with a view toward model selection and calibration. In particular, we 
consistently observe that the skew exhibits a power law of order between 0 and — 1 / 2 , which is in line with the skew- 
behavior of the models studied in Sections and while contradicting the behavior of purely continuous models, as 
well as model with a finite variation jump-component. Furthermore, we provide a simple calibration procedure for the 
index of jump activity of the process, T, which can be viewed as a new forward-looking tool to assess this fundamental 
parameter, complementing the popular rear-facing estimation methods based on high-frequency observations of the 
underlying asset’s returns. 


1.4. Outline. The rest of this paper is organized as follows. Section provides the probabilistic relationships on 
which we build our analysis, and introduces the class of tempered stable-like Levy processes. Section contains our 
results for the transition probability, volatility skew, and delta, under a pure-jump exponential Levy model. Section 
1^ presents the analogous results under a Levy jump model with stochastic volatility. Section contains numerical 
examples to assess the accuracy of the asymptotic expansions, as well as an empirical analysis of the short-term skew 
in recent S&P500 option prices. Section summarizes our results and draws some further conclusions. Finally, proofs 
of intermediary results are collected in the appendix. 
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2. Notation and auxiliary results 


Throughout, X := (Xt)t>o denotes a pure-jump tempered stable-like Levy process, as introduced [TB] and [Hj, 
defined on a filtered probability space (n, (^t)t>o,lP) satisfying the usual conditions. That is, X is a Levy process 

with triplet (0, &, p) relative to the truncation function (see Section 8 in [37]), where the Levy measure p is 

given by (1.6), for some constants C(1),C(—1) G [0, oo) such that (7(1) -I- (7(—1) > 0, T G (1;2), and a bounded 
function g : K\{0} —)■ [0, oo) such that g(x) —)■ 1 as a: —> 0. Let us also introduce the following additional technical 
conditions on the tempering function g, conveniently selected to facilitate the proofs of some of the results that follow: 

|g(a;) — 1 — < oo; (ii) limsup < oo; (hi) inf g(x) > 0, Vs > 0. (2.1) 

V |x| / |x|—i-cso 1^1 


(i) 


/|a:|<l 


|a:|<£ 


Here, a(l) and q;(— 1) are real-valued constants. We emphasize that the main results of Sections and only require 
condition (i) to be satisfied, which controls the behavior of the Levy density around the origin. In particular, a sufficient 
condition for all 1 < T < 2 is given by g(x) = 1-1- a{x/\x\)x + O(x^), as x —>■ 0. Next, dehne a measure transformation 
P —P, so that X has Levy triplet (0, b, v) under the measure P, where 


v{dx) = c{^^\x\ ^ ^dx, 

V|x|y 

is the T-stable Levy measure, and b is given by 


( 2 . 2 ) 


b:=b + 


'|x|<l 


!;(iy — iy)(dx) = b + C(l) / (1 — g(x))x ^ dx — C(—l) / (1 —(;(x))| 


I '^dx. 


In particular, the centered process {Zt)t>o, defined by 

:= Xt - jt, 

is a strictly T-stable process under P, and 


7 := E(Xi) = b + 


Y- 1 


-I-(7(1) J X (1 — g(x))dx — (7(—1) J \x\ {\ — g{x))dx. 

As is well known, necessary and sufficient conditions for St '■= to be a martingale are given by 

(i) / e^v{dx) < 00, (ii) b = — (e“ — 1 — xl|2;|<i)p(dx), 

and in that case, 7 can also be written as 

7 = - y {e^g{x) - g{x) - a;)C'(|y) \x\~^~^dx. 

By virtue of Theorem 33.1 in m, a necessary and sufficient condition for the measure transformation P 
well defined is given by 


(2.3) 


(2.4) 


(2.5) 


to be 


'X Mn 


( 2 . 6 ) 




where, hereafter, ip{x) := — Inq’(x). In what follows it will be useful to write the log-density process Ut := In as 


Ut = Ut+r]t:= f [ (p{x)N{ds,dx) +t f (e — I + ip{x))i>{dx), 
Jo Jmc, Jmc, 


which follows from Theorem 33.2 in m, and is valid provided that 

f — 1-I-y>(x)|£'(dx) < 00. 

J Mq 

We shall also make use of the following decomposition 


= 


xN{ds, dx) = 


xN{ds, dx) + 


0 JO 


xN{ds, dx) =: Z, 


(p) 




(2.7) 


( 2 . 8 ) 


(2.9) 


0 ^ —00 
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where, under P, and are strictly K-stable random variables with respective Levy measures 

v'^^\dx) := v‘''^\dx) ;= C(-l)|x|“^“4{2,<o}dT. (2.10) 

Finally, for future reference, we denote by Lz the infinitesimal generator of the process (.^t)t>o, which, for a function 
5 e C^, is given by 

{Lzg){x)=[ {g{u +x) - g{x) - ug{x))c(-^)\u\~^~^du. (2.11) 

Jmo ^ 1^1 ^ 

We conclude this se ction by collecting a couple of lemmas that will be needed in the sequel. The first one shows 
that the conditions in (2.1) are sufficient to justify the measure transformation P —>■ P, as well as the representation 

, and the proof therefore omitted (see further details in 


rip) 


(2.7). It is a simple adaptation of Lemma A.l in 


]): 


Lemma 2.1. Under (2.1), both (2.6) and (2.S) hold true. 


The next lemma provides key probabilistic relationships for the skew and the delta, that are fundamental to our 
approach. Concretely, let Rt := ln(5'(/5'o) be the log-return of the underlying, C(S'o,t, AT) = E (5'oe^* — Tf)"*" be the 
price of an option with strike price K, time-to-maturity t, and spot price Sq, and {So,t, K] a) be the price of 
the corresponding option under a Black-Scholes model with volatility cr. Then, the implied volatility a{K,,t) and delta 
A{K,t), parameterized in terms of the log-moneyness k := log{K/ Sq), are respectively defined so that 

dC{So,t,K) 


{S^, t, s^e^-d{K, t)) = C {So, t, Soe^), A{k, t) = 


dSo 


Jf=Soe" 


Finally, denote by $ and </> the standard Gaussian cumulative distribution and probability density functions. 
Lemma 2.2. Suppose that St admits a density function. Then, 

da{K, t) 


(i) 


e>P(S,>S„e-)-e-4.(-^iH)l^) 


dn 






(ii) A(A,t) = ^C(5o,L^oe") + e"P(.5t > ^oe"). (2.12) 

Jo 


Formula (2.12 -1) for the skew is well known in the literature (cf. [231 Ch. 5]) and is a simple consequence of the 


implicit function theorem, together with the identity 

dC{So,t,K) 

dK 


= -¥{St>K), 


(2.13) 


which holds true, e.g., under the assumption stated in the previous lemma. Formula (1.3) for the ATM skew is then 
obtained by using that a'{Q,t)'Jt —>■ 0 as t —)■ 0 (which holds under the framework considered herein [21], and, in fact, 
holds in a much more general model setting [3S|), and the s tandard approximations <l>(a;) = 1/2 + x/'/^ + 0{^) and 
1/{'\f^(j){x)) = l-\- x"^/2-\- Olx'^), as X —>■ 0. Formula (2.12 -ii) for the delta follows from the chain rule and (2.13), and 


formula ( |1.4[ ) for the ATM delta then immediately follows by taking k = 0. 

3. PuRE-juMP Levy model 

In this section, we study the short-time asymptotic behavior of the ATM implied volatility skew under the exponential 
Levy model St := Soe^*, where X := {Xt)t>o is a pure-jump tempered stable-like process as described in the previous 
section. For this model, a second order expansion for the ATM implied volatility is given in Theorem 3.1 of [2T], under 
a minimal integrability condition on q around the origin. Specifically, it is shown that 

d{t) = \/^aitY~2 y/^a2t^ + o{t2), t 0, (3-1) 


with (Ti := ¥,{Z((), where, under P, {Zt)t>o is a strictly stable process with Levy measure ^{dx) = C{x/\x\)\x 


-Y-l 


and 


0-2 := P(2'i < 0) C(l) / {e^q{x) — q{x) — x)x 


-Y-l 


dx 


dx. 


(3.2) 


— P(Zi > 0)C'(—1) / (e'^< 7 (x) — g(x) — x) |x| ^ ^dx. 


As explained in the introduction, the ATM skew is then related to the probability of the process X being positive. 
The following theorem gives an asymptotic expansion, in small time, for such a probability, which sometimes is termed 
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the positivity parameter of a process (cf. jHl P- 2 18]) . The expansion is explicit up to a term of order (see the 

I Iv J V (p) (n) J I 

subsequent Remark 3.5). Below, 7 is as in Eq. (2.4), ^ and are as in Eq. (2.9), and, finally, /^(p), /^(p), and 


fz, are the probability density functions of Z^^\Z{\ and Zi = z'f^ + Z^\ respectively. 

Theorem 3.1. Let X be a tempered stable-like Levy process with a Levy measure as described in Furthermore, 

assume that the condition 


l-Y-l 


dx < 00, 


/ - I - 

J\x\<i 

is satisfied for some constants a(l), «(—1) € R. Then, 

n 

P(Xt > 0) -P(Zi > 0) = +et^ +ft + o{t), t 0, 




where n max{A: > 3 : A:(l — 1/^) < 1}^ and 


\k-l 




k\ 


/ := 7(a(l) - a(-l))E(z(^)4(p) ( - 


\ < k < n, 

'' -L a(-l)E(Z.^”^ 


1 ^{z('’> + z<"’>0} 




+ E(-^i < O)C'(l) / {q{x) — 1 — a{Vjx)x 


-Y-l 


dx 


(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


— P(Zi > 0)C'(—1) / {q{x) — 1 — a{—Vjx)\x 


-y-i 


dx. 


Remark 3.2. The processes covered by Theorem 3.1 include stable processes, where q{x) = 1, and tempered stable 
processes as defined in [12], where q[x) = with a(l),Q;(—1) > 0. They are of particular 

importance for practical applications, and will be studied numerically in Section It is also important to note that 
condition (3.3) is the minimal condition needed for the expansion to be valid. That is, if (3.3) does not hold, the 
coefficient / is not well defined. 


If the Levy triplet (0, b, v) of X satisfies the martingale condition (2.5), then the previous result can be interpreted 
as an asymptotic expansion for ATM digital call option prices. Together with (1.3) and (3.1)-(3.2), it then gives an 
asymptotic expansion for the ATM implied volatility skew. Moreover, together with (1.4) and Theorem 3.1 of |2T|, it 
also gives an asymptotic expansion for the delta of ATM call options: 


Corollary 3.3. Let X be a tempered stable-like Levy process as in Theorem 3.1 with b and v satisfying (pi 
St := Sqc^* is a martingale. Then, 

(1) The ATM implied volatility skew satisfies 


da{K, t) 


dn 


27r /1 




P(Zl > 0) - ^ v) - (^e+ - (/+y)t + o(t)j, t^O, 




that 


(3.8) 


with tJi := E {Z^'j and 02 as in (3.2). 

(2) The delta of an ATM call option satisfies 


IL 

A(t) = V(^Zi > 0) + 'y ) ^'^-\-{(Ti-\-e)t'r-\-(^(j 2 -\-f)t-\- oit), t —^ 0. 


(3.9) 


k=l 

Remark 3.4. A few comments are in order: 

(a) It is important to point out that the leading order term of (3.4) is dit^~^^'^ for all Y G (1,2). Furthermore, the 
coefficients of (3.4) can be ranked as 




A ■ • • ^ 


Y^hft, t- 


0 , 
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where m = max{fc > 2 : k{l — 1/Y) < 1/Y} and, as usual, h{t) >- g{t) (resp. h{t) Y g{t)) if g{t) = o{h{t)) (resp. 
g{t) = 0{h(t))), as t —>■ 0. This can be compared to the expansion for ATM option prices given in Theorem 3.1 of [2T], 
where the first and second order terms are of order and t, respectively. 


(b) It is informative to note that the summation term in (3.4) comes from expanding the probability of a stable process 
with drift being positive. Specifically, we have 

n 

V{Zt + 7 t > 0) - P(Zi > 0) = ^ t ^ 0. 

k=l 

The other terms, e and /, arise as a result of the discrepancy between X and a stable process. In particular, this implies 
that the probability of X being positive at time t can be approximated, for small t, by the analogous probability for a 
stable process, up to an error term of order Similarly, the same probability can be approximated by that of 

a tempered stable process (as defined in Remark [T^ , up to an error term of order 0{t). 

(c) As mentioned in the introduction, the implied volatility of OTM options explodes as t —)■ 0, while for ATM options 

it converges to the volatility of the continuous component. Here we see that the ATM implied volatility slope also 
blows up as t > 0, with a sign that can easily be recovered from the model parameters. Indeed, when Z\ is symmetric 
(i.e. C'(l) = C'(—I)), it is of order with the same sign as the parameter 7 , i.e. the center of X under P, but, 

when C(l) > C(—1) (resp. C(l) < C'(—I)), it is of order with a negative (resp. positive) sign. 

Remark 3.5. There exist explicit expressions for E(Z(^) and f‘{Zi > 0) (see m and references therein): 


E (z+) = ^r(-y)Y 


(w) 


( 


cos ( — arctan ^ — tan 


'B 

.A 


Ytt 

~Y 


rii-y 




tan 


ttY 

) 




,11 (B /Yttw 

P(Zi > 0) = - + — arctan tan 


where A := (7(1) + (7(—1) and B := (7(1) — (7(—1). The derivatives ^^(0) can also be explicitly computed from 
the polynomial expansion for the stable density (see, e.g., Eq. (4.2.9) in [13]). Indeed, it follows that 




( 0 ) = (- 1 ) 


k— 


A 


—+1 
Y ^ ^ 

kiT 


where 


P = 


6 + Y 
2Y ’ 


6 = — arctan ^/3 tan f 


Ytt 

2 


sin(p/c7r)(^^) 

Co = cos ^ arctan tan 


/ttY 


and /3 = ((7(1) — (7(—1))/((7(1) + (7(—I)) and c = —r(—F) cos( 7 rF/ 2 )(( 7 (l) + (7(—1)) are the skewness and scale 
parameters of Zi. 


Proof of Theorem 13.11 

Step 1: Let A be a tempered stable-like process as in the statement of the theorem. In this step, we will show that 


(3.4) holds under the additional assumptions that the g-function of X satisfies (2.1 ii) and (2.1 iii), so that Lemma 2.1 


is valid. Throughout, we use the notation introduced in the previous section. Let us start by noting that 
P(At > 0) - ¥{Z, > 0) = + E((e-^‘ - l)l{z,>-^i}) =: Ii{t) + h{t), 

and we look at each of the two terms separately. For the first one, we have 

Ji(t)=P(Zi>-7ti-Y)-P(Zi>0)= f 

J — 

and, since fz is a smooth function (see e.g. [S^, Prop. 28.3), we can use its Maclaurin series expansion to show that 


^ fz{z)dz, 


(3.10) 


(3.11) 


N 




'(o)t“('-+) + J)). 




(3.12) 


For I 2 , we further decompose it as 

hit) = E((e-^‘ - l)l{z,>-^t}) + (e-''* - l)E((e-^‘ - l)l{z,>-7t}) + (e""* - l)E(l{^,>_^tj) 

=: I^it) + /|(t) -I- /|(t), 


(3.13) 
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(3.14) 


where it is clear that 

/2(t)=o(t), /|(t) = —77P (Zi > 0) t + o(t), t —>• 0. 

We use Fubini’s theorem on the first term to write 
I^it) = E((e — 1 + Ut)l{Zt>-jt}) ~ 

(e““ - l)P(^t > Ut < x)dx - / (e““ - l)P(Zt > -7^, Ut > x)dx - ¥.{Utl{Zt>-jt]) 

-00 J 0 

=: Ji(t) + J2(t) +J|(t). (3.15) 

Analogous arguments to those in m Eqs. (A.11)-(A.14)] can be used to apply the dominated convergence theorem 
and obtain: 

lim - Jj(t) = / (e“^ - 1 ) lim 7P(^t > Ut < x)dx = j {e~^ “ 1 ) / ^{<p{y)<x}’>{dy)dx =: r?!, (3.16) 


t —^0 t 


1 poo 1 poo poo 

lim = - I (e~^ - 1) lim jP(Zt > -jt, Ut > x)dx = - j {e~^ “ l) / ^{ip{y)>x}i'{dy)dx =: ■02- (3.17) 


t —>^0 t 


Finally, to deal with the third term of (3.15), we decompose Ut = fgf (p{x)N{ds,dx) as 


^t=J J {‘fiix)+a(^^^x)N{ds,dx) - J J a(^^^xN{ds,dx) =: U^^^ - U^'^\ 


so that 


J|(t) = +E([/f l{2,>-7t}) =: -JiHt) + 

First, for note that 


(3.18) 

(3.19) 


^^(p) 

\ J -jt^-Y -Zl’’’ 1 ' 

rip)\ 


— tvE(Zj + ytE(Z)^^/^(„)(—+ o{t), t —>■ 0, 


since sup/„(Ti)(z) < 00. Similarly, 


2^1 


=tvE(zri{2^>oj) +7tE(zr4(p)(-^r)) +o(t), t^O. 


From (3.20)-(3.21) and the fact that E(4”V2 (p)(“■^ i"'*)) = (“■^i^'*))> we get 




7 U) . 

rip) . 


rip)\ 


(3.20) 


(3.21) 


J|2(t) =4(a(l)E(4^)l{2,>0}) +a(-l)E(z4l{^,>0})) +7i(«(l)-a(-l))E(4^4 ^(„)(-Z 4 )) +o(t). (3.22) 
For J|i, we will show that 

=P‘{Ut^^i{Zt>-jt}) = dt + o{t), t —>■ 0, (3.23) 


where 


■d := C'(l)P(2'i < 0) / (a(l)a; — ln(7(a;))a; — C'(—l)P(Zi > 0) / (q;(— l)x — ln(7(x))|x| ^dx. (3.24) 

Jo J — 00 


Combining (3.13)-(3.24) then gives an asymptotic expansion for hit), which together with (3.10) and (3.12) yields 
(3.4), after some standard simplifications. To complete the proof we therefore only need to show (3.23). In order to 


do that, define /(x) := i^(x) + a(x/|x|)x, and, for e > 0, further decompose U^^ as 

41 )=/' j fix)N{ds,dx)= j f /(x)7V(ds,dx)+ / f /(x)lV(ds,dx) =: W, (3-25) 

Jo J Jo J\x\<.€ Jo J\x\>S 


Jl\t)=E{Ud^^\t)hz,>-^t}) +E{Ud^^\t)hz.>_r,t}) :=JiAt) + keit)- 


and let 


(3.26) 
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For future reference, recall that ip{x) = — lnq(a;), and 

J \f{x)\v{dx)< j |a^-j —-^x + 1 — q{x)\v{dx) + J |g(a;) — 1 — lng(a;)|u((ia;) < oo, (3.27) 

in light of (3.3), the boundedness of q, and the fact that (|2.1|) implies (2.8) as proved in Lemma 2.1 Now note that 
Ue ’ (t) is a compound Poisson process with drift; i.e., we can write 


N, 


(e) 




where := — f{x)v{dx), (-^t )t>o ^ counting process with intensity := J^^^^^i>{dx), and are 

i.i.d. random variables with probability measure We can also write 

T{e) ArCe) 


= 


atK^) 

P - f - ‘ f . 

/ / xN{ds,dx) = / / xN{ds,dx) +— t / xv[dx) =: 

Jq j Jo j\ x\<e „-_i J\x\>e „-_n 


|ai|>e 


and, under P, t v Z\ as t —>■ 0 (see 


, Proposition 1). Then, by conditioning on N\ , we have 
k,e(t) = = 0) + A(^)te-^‘'’‘E((/3(")^ + /(C^"^))l{ 2 ,>_^*}|iv(") = l) + o(t) 

= > -(^ + c(-))ti-^) + A(^)te-"'''*E(/(^(^V{4(=)+4(o>_(^+,<.,),}) +o(i) 

= d^^h + o(t), t^O, (3.28) 


where 


^(e) _ ^{e) > 0) ^ ^(s) E(/(e("))l 


{«r>op 

^ 1*00 ^ 

= C(l)P(Zi < 0) / (a(l)a; - lnq(a;)) x-^-^dx - ^(-l)? (Zi > 0) 


(a(—l)a; — Inq'(x)) |ai| 


and the second equality follows from standard simplifications. Moreover, in light of (3.27), 

1?—^0, e^O. 

For Ue^'^\t) we note that, by (3.27) and Theorem 10.15 in [21], we have 

\Ji,e{t)\ < < 2t / \f{x)\Hdx) =: ^ 0, 

J 

Finally, (3.23) follows since, by (3.26), ( |3.28 ) and ( |3.30 ), 

< liminf 


as e —>■ 0. 


(3.29) 

(3.30) 


t —^0 


< lim sup , 

t - t^o t- 


and the lower and upper bounds converge to i? as e —>■ 0 in view of (3.29) and (3.30). 

Step 2: Now assume that AT is a tempered stable-like process whose q-function satisfies (3.3), but not necessarily the 
additional conditions imposed in Step 1: (2.1-ii) and (2.1-iii). We would like to approximate it by a process whose 
^-function satisfies those conditions, and for which the result (3.4) is therefore known by Step 1. To do that, first note 
that since q{x) —>■ 1 as x —>■ 0, we can find eo > 0 such that inf|a;|<£(, q{x) > 0. Next, for each d > 0, let (n^^\ P*^"^^) 

be an extension of the original probability space (f2, P), carrying a Levy process independent of the original 

process X, with Levy triplet {0, given by 


{dx) := C”] o )e 


(A)® l{kl>Eo}|a;| ^ ^dx, := [ xv^l\dx). 


(3.31) 


In particular, R^^'> is a compound Poisson process and can be written as 




(a) 




(3.32) 
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where t>o is a Poisson process with intensity := ^R\dx), and are i.i.d. random variables 

with probability measure v^^\dx)/X^^\ Let us recall that, by the definition of a probability space extension (see 
[25), the law of X under remains unchanged. Also, all expected values in the sequel will be taken with respect 
to the extended probability measure P(‘^\ so for simplicity we denote the expectation under P^'^^ by E. Next, we 
approximate the law of the process X with that of the following process, again defined on the extended probability 


space (OW, J-(-5),p(-5)): 


:= X. 


i? 


(-5) 


(3.33) 


Then, the Levy triplet of is given by 

.= 5 + ^(<5)^ v'^^\dx)-.= c{^^\\x\~^~^q^^\x)dx-.= c{^^\\x\~^~^{q{x) + e~''^''/^l\^\>^g)dx, (3.34) 


so it is clear that satisfies the conditions in (2.1). Hence, the probability measure P can be defined as des crib ed in 
Section]^ using the jump measure of and note that for := E(X^'^^), we have, using the expression (2.4), 

^(5) = 5(5) + + c'(l) [ x-^{l-q^^\x))dx-C{-l) r \x\-^{l-q^^Hx))dx = j, 

y — ^ Jo J-i 

after plug ging in the above expressions for and q^^\ Now, since q^^'^ satisfies the conditions in Step 1, we know 
that (3.4) holds for X^^\ i.e. that 


lb 

P(xf^ > 0) -P(Zi > 0) = ^dkt^^^-^) +etv +/(‘5)t + o(t), t ^ 0, 


(3.35) 




where n := max{A: > 3 : A: (1 — 1/1^) < 1}, dk and e are independent of S and given by (3.5l-(3.6l, and 

:= 7(a(l) - z[^'>))+V{Z, < O)C'(l) / (q^^Hx) - 1 - ail)x)x-^-^dx 

1 Jo 

rO 

— P(Zi > 0)C'(—1) / {(f ^\x) — 1 — a{—l)x)\x\~'^~^dx. 

J —OO 


(3.36) 


Now, from the triangle inequality, it follows that 

P(xf ^ > 0) - |P(Xt > 0) - P(xf^ > 0)| < P(X( > 0) < P(xf^ > 0) + |P(X( > 0) - P(xf ^ > 0)|, (3.37) 

and, by conditioning on the number of jumps of the process we have 

:= |P(Xt > 0) - P(xf^ > 0)| = A('^)te^'''‘|P(Xt > 0) - P(X 4 + > 0)| + o(t), t 0, (3.38) 

which, in particular, implies that 

lim limf^TT.l'^^ = 0, (3.39) 


since e ^^^/^v{dx) —)■ 0, as (5 —> 0. By subtracting E {Zi > 0) + from the 


inequalities in ( |3.37 1, applying the expansion (3.35), dividing by t, taking the limit as t 
clear that 


0, and using (3.39), it is 


1 ‘‘ 

liin-(p(Xt > 0) -P(Zi > 0) -'^dkt^^^-^) - = lim/^'^). 

£->■0 t 5->-0 


Therefore, to conclude, it suffices to show that lim 

5->-0 

dominated convergence theorem. 


/, with / as in (3.7), which follows from (3.36) and the 

□ 
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4. Levy jump model with stochastic volatility 


In this section we consider the case when an independent continuous component is added to the pure-jump Levy 
process X. Concretely, let St := 5 * 06 ^*+^*, with X as in the previous section, while for the continuous component, V, 


we consider an independent stochastic volatility process of the form ( 1 . 10 ), defined on the same probability space as 


X, where [Wl)t>Q and are standard Brownian motions, relative to the filtration < P < 1; and 

a, 7 ,/!, and cr, are such that V and Y are well defined. Moreover, it is assumed that uq := cr(yo) > 0, and that there 
exists a bounded open interval /, containing i/Qj on which the function a is bounded, 7 and p. are Lipschitz continuous, 
and cr is a function. In the sequel, (j)s (resp. 0) denotes the probability density function of a A/’(0, 5"^) (resp. A/'(0,1)) 
random variable, while 


'I'(z) := f (j){x)d. 

Jo 


X. z G 


and ^ := 


4'(70 (*)' 


„1-Y 


dx. 


(4.1) 


Let us also recall that Lz, defined in ( 2 . 11 ), denotes the infinitesimal generator of the strictly stable process {Zt)t>o- 


The next theorem gives an asymptotic expansion for the probability of a tempered stable-like process being positive, 
in the presence of a continuous component satisfying the previously described conditions. 


Theorem 4.1. Let X be a tempered stable-like process as in Theorem \3. 1\ and V a diffusion process as described above. 
Then, 


V{Xt + Vt>G) = -+Y. +et^+ft'-^+ o(t^), 

k^l 

where n := max {A: > 3 : fc(l — Y/2) < (3 — T)/2}, and 

-L|4'(0), l<fc<n, 


t ^ 0, 


4 


k\ 


■= {lY hiVo) - |o-'(yo)7(?/o))</'^o(0)> 

^ ^ a(l)C(l) + a(-l)C(-l) _ ^ _ ^4(»)l(!,o)(l + y)))c 


Y - 1 

Remark 4.2. A few observations are in order: 


(4.2) 

(4.3) 

(4.4) 

(4.5) 


(a) Using the notation of Remark 3.4 the terms can be ordered with regard to their rate of convergence as 






t 


0 , 


where m := maxjfc : fc(l — Y/2) < 1/2}. A comparison of this and the expansion for ATM option prices given in 
Theorem 4.2 of [ 2 T], where the first and second order terms were of order and reveals that the convergence 


here is slower, as in the pure-jump case, unless (7(1) = (7(—1) (see (4.6) below). 


(b) As stated in the proof, there is another useful characterization of the dfe-coefficients in terms of a short-time 
expansion for a certain functional depending on ( |4.1[ ). Concretely, the coefficients di,... ,dn are such that 

E (^ (Zt)) = P{Zt + a{yo)Wl > O) - - = ^ d^t^ + 0(t"+i), t ^ 0. 

k=l 

When Zi is symmetric (i.e., when (7(1) = (7(—1)), it follows that E}!? (Zt )) = 0, and all the dfe’s vanish. In that case, 
the expansion simplifies to 


P (At -I- Ut > 0) = ^ + ft ’^ + o { t '^ 


t 0 . 


(4.6) 


(c) Interestingly enough, the correlation coefficient p appears in the expansion (4.2). Moreover, so does cr'{yo)j{yo), 
i.e. the volatility of volatility. This is in sharp contrast to the expansions for near-the-money option prices and implied 
volatility, given in Theorem 4.2 of |21j . where the impact of replacing the Brownian component by a stochastic volatility 
process was merely to replace the volatility of the Brownian component, a, by the spot volatility, a{yo). 

(d) Similarly, the leading order term di depends on the jump-component via Lz, and the parameters 7, 0 ( 1 ), and 
a(—1), also appear, containing information on the tempering function q, i.e. the Levy density away from the origin. 
This is again in contrast to what was observed in Theorem 4.2 of |21j . where the leading order term only incorporated 
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information on the spot volatility, cr(yo), and the approximation was altogether independent of the g-function. In short, 
the ATM skew is more sensitive to various model parameters, which may be anticipated since it is a measure of the 
asymmetry in the volatility smile, i.e. the difference between the volatilities of OTM call and put options, while the 
ATM volatility is a measure of the overall level of volatility. 

(e) Tempered stable-like processes are a natural extension of stable Levy processes. In the pure-ju mp case, the 
“deviation” of X from a stable process does not appear in terms of order lower t han (see Remark 3.4-b). Here, 
recalling that Xt has the stable representation := Zt + yt under P, Theorem 4.1 implies that, 

a(l)C(I) + a(-I)C(-I) 


' (At + Hi > 0) - P (Af “ + H > 0) = 


Y - 1 


pOO ^ ^ ^ ^ 

■ / cj)„g{x)x^~^ dxt^^ t 

Jo 


0 . 


In other words, for small t, one can explicitly approximate the positivity probability of Xt + Vt by that of Aj**’* -|- H, 
up to a term of order higher than 

(f) One can find a more explicit expression for the constant / by noting that 




4^(70 {x). 


i<^iyo)) 


l-Y 




-EIVLi 


(c^( 2 /o)) 


l-Yn-Y+l 
^ ^ 2 2 




using the well-known moment formula for centered Gaussian random variables. Moreover, we can further show (see 
[55] for the details) that the first two coefficients in (4.3) are given by 


:= 




{a{yo)VY Jo 
1C'2(I)-C'2(-I) 

2 {aiyoWyY^ Jo Jo 


pOO 

/ {J>{x)-m)^~^dx = -{C{l)-C{-l)) 

Jo 


{a{yo))-^2-Y 




^oo /*oo 


Y^YiY - 1) 

{(x + y)(t>{,x + y)- x(j){x) - #(i/)) {xy)~^dxdy. 


(4.7) 

(4.8) 


In the case when the asset price St = Soe^*-~^'^* is a well defined P-martingale, the previous result can be viewed as 
an asymptotic expansion for ATM digital call prices. Moreover, for models of this form, Theorem 4.2 in [21] supplies 
the following short-term expansion for the ATM implied volatility, 


where 


a{t) = a{yo) + (Tit ^ +o{t ^ t 

(C(l) + C(- 1 )) 2 -^ 


0 . 


CTi := 


y(r-i) 


r ( 1 - ^ ) (T{yo) 


l-Y 


(4.9) 


(4.10) 


The above results can then be combined with (1.3) to obtain an asymptotic expansion for the ATM implied volatility 
skew. Additionally, using (1.4) and the ATM option price expansion in Theorem 4.2 of m , we also obtain an asymptotic 
expansion for the delta of ATM call options: 


Corollary 4.3. Let X and V be as in Theorem 4-1 with b and v satisfying the martingale condition (2.5), y, = 
and a, 7 , and a such that (e'^‘)t>o is a true martingale. Then, 

(1) The ATM implied volatility skew satisfies 


9(t(k, t) 


Ok 


K — 0 


IL 

VJjT dk tO 




fc=i 


-1- I V 27r 1 

(To 


ldi]t 


i-i 


-|- o{t 2 )^ f — 


where c := 7 — \p(t'{ yo)l{yo)j (rnd, if C{1) = C(—1), the expansion becomes 


da{K,t) 


Ok 


K=0 o'o 

(2) The delta of an ATM call option satisfies 


= - \-(y/lnf +-aA t 2 -(-o(i 2 )^ t —)• 0. 


A«) = 5 


k—1 ^ ^ 


(4.11) 


(4.12) 


(4.13) 
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Remark 4.4. The previous result shows that the order of convergence and the sign of the ATM implied volatility 
slope can easily be recovered from the model parameters. In the asymmetric case, i.e. when C{1) ^ C'(—1), it blo ws up 
like and has the same sign as (7(1) — (7(—1). However, when (7(1) = (7(—1), the summation term in (4.11) 

vanishes, and the slope converges to a nonzero value, — c/cto, as t —>■ 0, as in jump-diffusion models. In both cases it 
is observed that for a fixed value of the index of jump activity, Y, the short-term slope is less explosive than in the 
pure-jump case. 

Proof of Theorem 14.1[ 


Step 1: We first show that (4.2) is true when the functions /i and a are assumed to be bounded, by considering the 
stopped processes 

:= ^(Fmt), dt := cr(rtAr), r := inf {t : T* ^/} , (4.14) 

where / is a bounded open interval containing such that a is C^, y, and 7 are Lipschitz, and a is bounded on I. 
Note that due to the continuity of y and <j around i/Oj we can find constants 0 < m < M < 00 such that \yt\ < M and 
m < at < M. Throughout the proof, we set uo = a'{yo), yo = y{yo), <to = C({yo), 70 = 7(1/0)) and p := \J\ — p^. As in 
the pure-jump case, we also start by assuming that the g-function of A satisfies ( 2 . 1 -ii) and ( |2.1| -iii). Then, the idea is 
to reduce the problem to the case where yt and at are deterministic, by conditioning the positivity probability on the 
realization of the process (lT(^)o<t<i- To do that we follow similar steps as in the proof of Theorem 4.2 in [21]. On 
a filtered probability space (0,.F, (.Ft)t>o,F) satisfying the usual conditions, we define independent processes X and 
such that the law of (At)o<t<i under P is the same as the law of (At)o<t<i under P, and (lT(^)o<t<i is a standard 
Brownian motion. Also, for any deterministic functions y := (/1s)sg[o,i]) ^ ■= (^s);i 6 [o,i]) and q := (gs);sG[o,i]) belonging 
to (7([0, 1]), the set of continuous function on (0, 1), we define the process (W’‘^’'^)o<t<i as follows: 




yudu + pqt+ p / 
Jo 


0<t<l. (4.15) 

With this notation at hand, we consider a functional $ : [0,1] x (7([0, Ij) x (7([0,1]) x (7([0,1]) —)■ [0,1], defined as 

^{t,y,d,q) -ViXt + Vr’^^O). (4.16) 

Then, for any t S [0,1], 

P( 7ft + Vt > 0| Wl , s G [0,1]) = 4)(t, (/ls)sG[o,i]) (cts)sg[o,i]) (9s)sg[o,i])) (4-17) 

where qs '■= Jq dudW^. For simplicity, we omit the superscripts in the process unless explicitly needed. 

Throughout the proof, we assume that y and d satisfy the same uniform boundedness conditions as y and it; namely, 
m < dt < M and \yt\ < M for any t G [0,1]. 

As in the pure-jump case, we define a probability measure P' on (17, A), analogous to the probability measure P 
described in Section]^ but replacing the jump measure N of the process X by the jump measure of X. We also define 
the strictly stable process Zt := Xt — jt, where 7 := E'(Ai), and E' denotes the expectation with respect to the 
probability measure P'. Note that the law of {Vt)t<i under P' remains unchanged and, under both P and P', 

r^Vt ^N{t^yl +r^pqt,p^{d*t f) (4.18) 

where, for t G (0,1], 


^^t ■= 


[ ysds G 
Jo 


a. := 



dgds G [m, M], 


(4.19) 


Now, note that P (Aj -|- V) > 0) = P(t ^Xt + t ^Vt > O) converges to 1/2, as t —)■ 0, by Slutsky’s theorem, and the 
facts that t~^l'^Vt A ^ A/’(0,CTq) (cf. [2T], Eq. (4.47)), and Xt Z\, where Z\ is a strictly F-stable random 


variable (cf. 


], Prop. 1). In order to find higher order terms in the expansion, we investigate the limit of the process 

1 


i?t:=P(At + yt>0)- 


2 ’ 


(4.20) 


as t —>■ 0. In terms of the functional 4), Rt can be expressed as 

Rt = E^P (At -I- Vt > 0| wl , s G [0,1]) ~ 2 ) ~ (Ms)sg[o,i]) (d’s)sG[o,i]) (9 s)sg[o.i]))) 


(4.21) 
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where $ {t, ji, a,q) = ^ {t, fl, a, q) — 1/2. Note that 

= E'fl 


{t 2Vt>-t 2Zt-jt2} 


- 1 , 


{t 2Rj>-t 2Zt-jt2} ^{VFl>0}, 

=: fl, a, q) + hit, fi, a, q). 


+ E'((e-^‘-l)l^,_ 


{t 2Vt>-t 2Zt-^t2}l 


(4.22) 


and we proceed to analyze the two terms separately. For the first one, we use (4.181 to show that 
Ii{t,fl,a,q) =E'( 


1 _i (t>pat{x)dx) 
-t 2 Zt-(^ + ill)t2 -t 2pgj Z 


= E' 


■E' 


-t 2 Zt-{^ + po)t2 -t 2 pq^ 

-t~2 Zt-(^+il^)t2 -t~2 pqt 
-t-2 Zt-t-2 pqt 


(l)p!i*ix)dx) +E' 


-t 2Zt-t 2pqt 

_ 1 „ , .1 _i {(t>p)y*ix) - (fpaoix)) dx) 

-t 2 Zt-[j+pojt2 -t 2 pgt 


^ f>paoix)dx)+E' 

I-t 2 Zt-(^+po)t2 -t 2 pqt 

=: llit,fl,a,q) + I^it,d,q) + lfit,q) + I^it,d,q), 


-t 2Zt-t 2pqt 


4>pP*{x)dx) 


(4.23) 


where we recall that denotes the density of a N{0,(t‘^) random variable. For the first term of (4.23), we have 

E |j/(t,/2,d-,g)| < t^p„(0)E|/2/ -/Jol < t^pm(O)- / E\p,s - pio\ds = 0{t), t ^ 0, (4.24) 

^ Jo 

where the last step follows from Lemma [A.l (i) below. Similarly, by Lemma A.l (ii) and the easily verifiable fact 
sup If^pafix) - (t>pao{x)\ = |/><T*( 0 ) “ (faoiO)^ 




E|/i(t,cr,(?)| <t2|7 + Mo|E|/)pff*(0)-(/)p„o(0)| < 


^H7 + MoI 

'■pV^ 


m- 


E |crt - CTol = 0(t), 


(4.25) 


For the third term I fit, q), assume 7 ' := 7 + /io > 0 (the analysis when 7 ' < 0 is identical), and write 

r-t~2Zt 


t 2 /i(Lg) - y</^o( 0 ) = ^ "E' 


+ t 2 E' 


-t 2Zt-^'t2 
-t~2Zt 


^{(fpaoix-t 2 pa^wO - (j)paoi-t 2paQwO)da 


^{4>paoix-t 2 pq^) - - t 2 paQwO)dx\ 

l-t 2Zt-l't2 ' 

+ y(</>pao(-i"®PCroW^) - f/aoW) 

=: /i’^(t,«;^) +/f’^(t,(7,w;^) +/f’^(<,'ui^), 

where ish € K. First, for w^), use Fubini’s theorem to write 

/ oo 1 ^ 

{4>paoix - paot~2w0 - (fpaoi-P<Xot~ ^ wO) Mx)dx, 

-00 


(4.26) 


where Jtix) := P'(—ts ^ x — 7 '!^ ^ < —t 2 y x) < kt 2 


\-Y-l 


with the inequality being a special case of 


(4.35) below, and holds for some constant k, all a; ^ 0 and t < 1. Moreover, Jt{x) ^ yC'(—as 
t —> 0, which follows from fzix) ~ C(a;/|a;|)|a:|“^“^, \x\ —)■ 00 (cf. [37], 14.37). On the other hand, it is easy to see 
that E (^<j)pryg{x — paoWf) — c/p^o (—pcrobFi)) = 4’<yo{x) — </cro(0) = 0(a:^), as cc —>■ 0, so, in light of the above relations, 
the dominated convergence theorem can be applied to E(/j’^(t, IF/)), to obtain 


poo 

lirnt-^E(//’i(t,VF/)) =f (C(l) + C(-1)) / (/).„(cr) - </.„(0)) 

Jo 

For the second part of (4.26), we can find a constant k such that 

1 / / / _i '-'1x1 ~\ —I ^ _1 <j\ 

\(Paopix-t v*?*) - 0crop(a; - f ^pcroW)\<K\t 2 q^ - f 2aow\, 


(4.27) 
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bb/)) — 0(f2) _ o(t 2 f —>• 0. 


(4.28) 


Finally, = 0 since E ((/ip^o (—p^obFi)) = (f’croiO), which together with (4.26 )-(4.28) gives 

poo 

=7'<?i’rTo(0)t5+y (C'(l) + C'(-1)) / {(j)ao{x) - (j)ao{0))x~^~'^dxt^ + o{t^), t ^ 0. (4.29) 

Before handling the fourth part of (4.23), let us introduce some further notation. First, note that by Ito’s formula, 

qt= [ dsdW^ = aoWl + f f d'^judW^dWl + f f (K^u + ^d'^^DdudW^ =: aoW^ + (4.30) 

Jo Jo Jo Jo Jo ^ 

where = a'(F„)l{„<.r}: = ct"(Fu)1{„<^}, := a{Yu^T)^ and % := 7(F„Ar)- Also, define 

CJ ■■= a'^i^dW^^dWl = id'7o {{Wlf -t)^ {{Wlf - l) 

and, for reals and let 
lt{i^d,q) = E' 


(4.31) 


-t 2 Zt-t 2 p(^ao-w^+i^ , ^ ^ r-t 2 Zt-t 2 paQil)^ 

(t)pa-{x)dx) +E'[ ^ (l)ps*{x)dx) 

' ^ J — t 2 — f 2 n( (Tn-j/'a -Uf I ' 


-t 2 Zt-t 2 pqt 
0 


. E' 


— t 2 Zt — t 2 pao 


-t 2 Zt-t 2p(ctoUi1+{) 

1 (l)paoix)dx) +E'( {4)p^*^{x) - (i3p^^{x))dx\ 

2 P<7qW^ ' ^ J — t 2 Zf — t 2 p(TQW^ ' 


= J^{t,d,q,w^,^) + J^{t,d,w^,^) + + J'^{t,d,w^). 


For the first term, by Lemma A.l (v), 

E\j^{t,d,q,Wt\Cl’°)\<MO)pt~H^^t\+^^l-JJ\)=Oit), t^O, 


(4.32) 


(4.33) 


For the second term, Cauchy’s inequality, (4.31), and Lemma A.l (iii) can be used to show that E|^y (^pg.*(0) — 
'/’pffo(0))| = 0{ti), as t —>■ 0, and, thus, using again that |(/)pg*(a;) — ^pCTg(x)| attains its maximum at x = 0, 


E(j2(t,d,lF/,et’°)) =e(e'( J 


— t 2 Zf— pCToW 


— t 2 Zt — po-ow — t2 ^a-g'Yo{w^ — l) 
Z) 


(j)p^„{x)dx\ )+0{t), t^O, 

/ w—Wl/ 


where above we also used that {^l’^ ,Wl) = (^crgyot ((bF^ )^ — l) , )• Thus, by the dominated convergence 

theorem, 


limt-5E(j2(t,a,bF/,e^°)) = ^-a'^l^E{{{Wlf - l)cj,p,,{pa^Wl)) = -'^</>.o(0)a'7o, 
where the last step is the result of elementary calculations. To find the second order term of J^, define 

_ i - 

P — t 2 Zt — pCTQW 
I -t~ 2 Zt-pcrow-t2 f(To7o 




(4.34) 


J^(t, w) :=t 2 E' 


4>pa-oix)dx'^ 


/ OO 

{(t^aopix + pCfQw) - 4>cop {pCfQw)) Jt{x, w)dx, 

-OO 


- - l) (t>pao (pcroih) 


where 


Jt{x,w) :=¥'{t "^Zt<x<t 2 Zt + a ^crp 7 o - 1)) - P'(t " |crp 7 o - l) < x < t Zt) 




i _ i 

rt2 Y X 


= /1 i , i „ fz{z)dz, 

Jt 2 Y x — l Y ^( 7070 ( 111 ^ — 1 ) 

for which, we can use the estimate /z(z) < R\z\~^~^ (see, e.g. [37], (14.37)) to find a function / such that E (/(IF^)) < 
OO and, for all x ^ 0 and < < 1, 


\Jt{x,w)\ < f{w)t 2 |x 


-Y-1. 


(4.35) 
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see [33] for the details. Next, by applying the dominated convergence theorem twice, we get 

/ OO 

E ((<^o-op(a^ + pcToVbi) - (l)aop (pcToVbi)) Jt{x, W^)) dx 

-OO t^O / 

= f <^o7o J E {{(j)aap{x + paoWl) - {pcjqWI)) {{Wlf - l)) C'( j^) \x\~^~^dx 

r? r°° .7.2 

= ((7(1)+ C'(-l))^cro7o / {(j)ao{x){^- 1 )+(l>ao{Q))x-^-'^dx, 

^ Jo ^0 


l-Y-1 


(4.36) 

(cf. |37|, 14.37), and 

standard calculations. Note the second application of the dominated convergence theorem above follows from (4.35), 
and the boundedness of ^ipo-o- The first application of it can also be justified using (4.35) for |a;| > 1, but for \x\ < 1 
we use Taylor’s theorem to switch the order of limit and integration. More precisely, we can write 


where the final two equalities follow from the tail estimate fz{x) ^ C{x/\x\)\x 


1 


(l)aop{.x + paoib) = (t)aop{paQw) + x (t>'c„p{x + pcrow)L n + x a; Cop(* + U=e 


(4.37) 


where 0 < |^a;| < |a;|, so 


E 


( J ^ {{<t>aop{x + paoWl) - (t)aoP (pcTobbi )) Jt{x, W^)) dx'j = J ^ E(^i X^<j 3 ”^p{x + pcroWl)\^^^^ Jt{x, WD^dx, 


becaus e (j)'„^ p{x + paow)\^^^ = -{pw)/{pao)<paop{po'ow) and E {Wl(j)cropipaoWl)Jt{x,Wl)) = 0, due to symmetry. 

Then, (4.35), and the fact that (j>" is a bounded function, allows us to apply the dominated convergence theorem, 

limf-V^ J ^E(^y Cop(a; + pcrou;)|^^^^ Jtix,Wl)^dx = J ^ limt-^E( (<))<,„p(a; + pcro^Tj^) - <(><top (pcroW^)) Jtix,Wl))dx, 

where we have again used (4.37). Finally, from ( |4.34 ) and (4.36) we get 

0 / / \ \ poo 

^^^E (J2 )=-</>.„(0)t5-(C(l) + C(-l)) / ((/..„(x)(0))x-^-ida:t^ 

P^<Xolo \ \ J J Jo 

1 3 Y 3 Y 

+ ((7(1) + (7(—1 ))—2 / (j)cro{x)x~'^~^^dxt^~ +o(t^~), t—>■ 0. (4.38) 

'Xq Jo 

Next, for the third term in (4.32), let P denote the probability measure on (7't)t>o) defined in Sectionj^and let 

Z := (Zt)t>o be the process defined in (2.9). Note that by the independence of Z and W^, and the fact that the law 


of Z under P' is the same as that of Z under ) 


E(J3 =E 


-t 2Zt 


= E 




<(’ (y) dy^, 


(4.39) 


where the last equality we used the self-similarity relationship = Zsf Therefore, it is sufficient to find the 

asymptotic behavior, as t —>■ 0, of E(4'(Z4)), with 4' as in (4.1). But, since 'P(z) has continuous and bounded 
derivatives of all orders, an iterated Dynkin-type formula (see mr Eq. (1.6)) can be applied to obtain 

^ ^ j-k +‘71+1 /*1 ^ 

E (^(Zt)) = ^(0) + V ttT|^( 0) + ^ / (1 - a)”E (L’^+i^(Z„t)) da, 

, n'. Jo 


k=l 


for any n G N, where Lz is the infinitesimal generator of the strictly stable process Z, defined in (2.11). Therefore, 

“ Uo-^^£|vl/(0) V 


E(+ {t,W}))=Y^ 


k\ 


(n+i)(2-y) , 

'’+0{t 2 y t->■ 0. 


(4.40) 


Finally, for the fourth part of (4.32), use Taylor’s theorem to write 

+a*(2^) = <^P<To(x) - _ cTp) h:^(dt)(dt - cro)^ 




p^ao' 
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where —)■ 0, as > cto, and the boundedness of CTj away from 0 and oo allows us to find a constant K such 

that 0 < \hx{a'^)\ < K, for alH < 1 and all a; S K. From the latter, and Lemma A.l parts (hi) -(iv), it follows that 

c — t Zt 


E 


{At,a,W,^))=-E[{at-ao)E'[[ ^ ‘ ^ (i _ I ) + 0(t) 

V ^ J-t ~2 Zt-t ~2 paow^ ^0 P ^0 ^ / 






—t 2pCQl 


— t 2 Zf — t 2 p(ToW^ 




where the second equality follows from Lemma A.l (vi). To handle the last expression, let us first note that, condi¬ 
tionally on Wf, Jq Wlds is normally distributed with mean and variance, tW}/2 and t^/12, respectively. Therefore, 
again using the probability measure P and the process Z := (^t)t>o, as in (4.39), we can write 

.{J^{t,d,wl)) = e{U--^wI 


-t- 2 ^Ei 


c^o7o 


r-t 2Zt 
' —t~2 Zt — p(7ot~ t Wl 


i - ^)dx) = ^E(i>{Z{P ^))), 


where we have used the self-similarity relationships s 2 = Wl^ and s^Zt^Z^t, and the notation 


^( 0 ) := e(vFi^ 


fZ+paoWx^ a;^ \ 1 

<(>P^o(2^)(l - T^)da;) = —(l)„„{z){al - z^)p‘^p. 

z P CTo ^ ^0 


Since ^( 2 ) has continuous and bounded derivatives of all orders, we proceed as in (4.40) and obtain 

-2-^E (J'^ d,q, Wl)) = iL(0)t5 -I- Lz4'(0)t^ + o{t^), t -)■ 0, 
o-o7o 

where ^( 0 ) = pp^4>(fi) and 


(4.41) 


r °° 1 

Lz'^{Q) = ((7(1) + C{-l))pp^ / { — - V?) - cro«!><To(0))R"^"^dw, 

Jo o-Q 

which trivially follows from (2.11). Thus, combining ( |4.32 )-(4.33), and (4.38 )-( 4.41| ), gives an asymptotic expansion 
for E (/(‘(t, (T, g)), which, together with (4.23l-(4.25l and ( |4.29[ ), finally gives 

E(/i(t,/2,d,g)) = -p (f - |cro7o)((>„o(0)t5 

k^l 

C(l) + C(- 1 )_ p. 


alY 


poo 

( 7 '- |cro 7 o (1-f y)) (()<,o(x)x"’^+^da;t^-f o(t^), (4.42) 


as t —>■ 0, where n := maxjA: > 3 : fc(l — Y/2) < (3 — T)/2}, and we have used integration by parts to write 

poo -I poo 

/ {(l^aoix) - (j)aoi0))x~'^~^dx = -^ (x)x"’^+^(ix. 

Jo ^ 0 ^ Jo 


Now consider the second part of (4.22). By using (2.7) and similar steps as in (3.15), followed by the decomposition 
(3.18), we can write 

l 2 (t, u.d.q) =—E', 1 - 1 - 1 ,)-|- 1 - 1 - 1 ,)-|-i?(t, u, (f, g) 

\ * {t~2Vt>-t~2Zt-jt2}J V t {t~2Vt>-t~2Zt-jt2}J 

= -ii(t, p., d, q) + hit, p, d, q) + R{t, p, d, q) (4.43) 

where E{R{t,p,,d,q)) = 0{t), as t —?> 0, and E{ii(t,p,d,q)) = 0(t) since is a finite variation process. We further 
decompose hit, fl, d, q) as 

T't 7(^)1 


=: a(l)l 2 {t,fl,d,q) + a{-l)l 2 {t,p,,d,q). 


hit,u,d,q) = ail)E' iz}^ 1 ^ 1- i/'-fpi 1 ,) + Q!(—l)E^(Zi'" 1 ^ 1- i/'-fni “(ni\ 1,) 

V ; V t {(-2 V^>_t-2 (Zj^'-t-^i M-7t5}^ ' \ J \ t {(-2 Vj>_t-2 (z)P>+z)"M-7i2}/ 


(4.44) 
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and we look at the two terms separately. For the first one we have 




= E' Z, 




'(p) 


-t 5 ipijt 

-t 2Zl^^-t 2pg, 


-t 2 2 pq-t J-t 2Z),^>-t 2 pqt J-t 2 pqt 


1 4>pGt{x)dx 


/ -t 2 pq^ „oo 

+/ _i ](l)pat{x)dx\ (4.45) 

-t 2 Z:^ —t 2/504 J —t 2 OQt ' ' 


+ 


and the third integral is zero since = 0. For the first one, use the independence of Z'^'^ and to write 


(p) 




E' 


Kp) 


r-t-2Z[^^-t-2pq, 


'-t 2 (z)P4z)">)-(7 + Pj)t2-t 2 


^ (t>pai{x)dx <t^E'|z(P^|(^p^(0)(t^-5E'|z(")| + (|7|+M)t5) 


= 0{t^ •2) = o{t ^ 0. 


(4.46) 


Finally, for the second integral, let iu^ € M, and write 

nt~ 2 pcTQW 




t 5 2)”'+* 2pqt 


t 2 pqt 


^pij^{x)dx^ =E'[z'f'' ^pif^{x)dx^ +E'(^Z^^ 

J t 2 pq^ 

~ / w/ \ /■* 2Z[^^+t 2 pqt 

+ E'(zlP^ , (l)pa-ix)dxj 

^ Jt 2.zlP>+t 2 panw^ ' 


t 2 Z^ -\-t 2 pcTQW^ 


(j)pa’;{x)dx 


ft 2 pcToW^ 


ft 2 Z^^ -\-t 2 paoW^ 

=: Ji{t, d, q, w^) + J 2 {t, d, w'^) + Jsit, d, q, w^), 

and observe that 

E\Ji{t,d,q,Wt) + J3(t,d,q,Wl)\ <2p4>pmiO)t^~^E' z[^'’ E|gt - crobF/1 = 0(t), t ^ 0, 


(4.47) 

(4.48) 


by Lemma A.l (v), which implies that E \qt — croVF/| = 0(t^/^). Next, write 

pOO ^ 

J2{t,d,w^) = W paow^)dy 

- j_^ -i z[-^<y<3})^P-t (y - t~h<xow^)dy 

=: J2{t,d,w^) - J2{t,d,w^). 

It then follows that 

limt"^E(j](t,CT,IF/)) = y^~^(l>aq{x)dx, E (j|(f, d, VF/)) = o(t^), t ^ 0. 


(4.49) 


(4.50) 


The first relation above follows from a similar procedure as in the proof of (4.26) in [21], while the second holds because 
the jump support of Zl ' is concentrated on the positive axis. Combining (4.45l-(4.50l gives 

m , 


E {I^{t,p,a,p)) = 


^l-Y 


Y -1 


4>cro{x)dxt 2 +o(t 2 )^ t —>• 0. 


Finally, for the term /|(f) in (|4.44[), the same procedure can be used to obtain 


E (/|(t,^,cr,p)) = j ^</)a„(x)dxt''2"' +o(t"2'^)^ t^O, 


which, together with ( 4.43| )-( [4(44 ), yields 

a(l)C(l)+ a(-l)C(-l) , T_v , ^ 

- / (^,jg(x)x^ ^dxt 2 -(-o(t 2 )^ f —0. 

Jo 


E{l2{t,p,a,q)) = 


Y -1 


(4.51) 


(4.52) 


(4.53) 


Combining ( |4.20[ )-( |4.22[ ), ( |4.42[ ), and ( |4.53[ ), then gives ( |4.2[ ). 

Step 2: The next step is to show that the expansion (4.2) extends to the case when the g-function of X does not 
necessarily satisfy conditions (2.1ii) and (2.1iii). That can be done exactly as in Step 2 of the pure-jump case, by 
defining a process X on an extended probability space (fl,.F, P), satisfying those conditions, and using the triangle 
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inequality to show that the terms of order lower than t in the asymptotic expansion of P(Xt + Vt > 0) extend to 
P(Xt + Ft>0). 

Step 3: Lastly, we will show that the expansion extend s to t he general case when a and fi are not necessarily 
bounded functions. To that end, define a process (Vt)i<i as in (1.10)), but replacing cr(Yt) and with the stopped 
processes at := a{Yt/\r) and fit ■= f^{Yt/\T), introduced in (4.14). By Steps 1-2 above, the asymptotic expansion (4.2) 
holds for the process X + V. For it to extend to the process X + V, it is then sufficient to show that 

|P {Xt + -Ft > 0) - P (Xt + Ft > 0) I = 0(t), t ^ 0, 

because (3 — y)/2 < 1 for Y G (1, 2). But since Vt = V) for t < t, we have 

\¥{Xt + Vt >0)-¥{Xt + Vt>0)\ = \FiXt + Vt>0,T <t)-F{Xt + Vt>0,T <t)\< 2P(r < t) = 0(t), 

as t —>■ 0, where the last step follows form Lemma 4.1 in [2T]. □ 

Remark 4.5. We hnish this section by briefly considering the OTM implied volatility skew, under an exponential 
Levy model St := , where (W)t>o is a Levy process with generating triplet (0,6, u), and (Wt)t>o is an 

independent standard Brownian motion. To this end, fix a log-moneyness level k ^ 0, and recall the following well 
known asymptotic relationship (see, e.g., [37], Corollary 8.9), 


P {Xt -b aWt > k) = tv{[K, oo))l{^>o} + (1 - tu((-oo, k]))1{„<o} + o{t), t 0, 
and that for the implied volatility of an option with log-moneyness k. Theorem 2.3 in m states that 


(T^(A,t)t — ^1-b Vi(t, At)-b o(^-^)), t —>■ 0, 


2 In ^ 

under some very mild conditions on the Levy measure, where 


T, ^ 1 1 ,^4v^ao(At)e 




(4.54) 


(4.55) 


(4.56) 


and ao(At) := — e”)+u((ia;)l{K>o} + ~ Plugging the above relations into (2.12 4) then 

yields the following higher order expansion for the OTM implied volatility skew: 


da(K, t) 
Ok 


1 At 

2t In - = — 

t At 


Vi(t, At)\ At 1 




lAtl 2Ini 


V 1 




t —y 0, 


(4.57) 


where 6o(At) := —e'‘(u([At, oo))1{„>q} -b u((— oo , At])l{K<o}) (see Chapter 3 in [35] for the details). It is noteworthy that 
this expression can be obtained in much more generality than the ATM skew expressions in Corollaries |3.3| and |4.3| 
For instance, the above result does not depend on the presence of a nonzero Brownian component, or whether X has 
hnite or infinite jump activity. Furthermore, it is interesting that this expression for the limiting skew can also be 
deduced (at least formally) by differentiating the expression (4.55) for the implied volatility - something that is not 
possible for the ATM skew. 


5. Numerical examples 

In the first part of this section we use Monte Carlo simulation to assess the accuracy of the asymptotic expansions 
presented in Sections [^ and [^ In the second part we investigate the implied volatility skew in S&P500 option prices 
and how it compares to the short-term skew of the models in Sections [^ and [^ 

5.1. Accuracy of the asymptotic expansions. In this section we carry out a numerical analysis for the popular 
class of the tempered stable Levy processes, with and without an independent continuous component (see, e.g., [5], 
[12], [27]). They are an extension of the CGMY model of [lO], and characterized by a Levy measure of the form 

p(dx) = 


(5.1) 
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where C'(l) and C'(—1) are nonnegative such that (7(1) + (7(—1) > 0, (7 and M are strictly positive constants, and 
Y G (1,2). The martingale condition (2.5) also implies that M > 1. Note that in terms of the notation of Section]^ 
we have 0 ( 1 ) = —M and a(—1) = (7, and the constants 7 and rj defined in Eqs. (2.4) and \2.7\ are given by 

7 = -r(-y)(C(l)((M - 1)^ - M^) + C(-1)((G + 1)^ - G^)), 

r; = r(-r)(G(l)M^ + G(-l)G^). 

In the pure-jump case. Theorem |3.1| and Corollary |3.3| present short-term expansions for digital call option prices and 


the implied volatility skew. For tempered stable processes, the dfc-coefficients are as in (3.5), while 

e = (5.2) 

/ = -7(m + g)e(4^V2w + r(-y)(P(Zi < o)g(i)m^ - p(Zi > o)g(-i)g'"). (5.3) 

It is informative to note that the terms can be further simplified in the CGMY-case, i.e. when G := G(l) = G(—1). In 


that case, which implies > 0}) = E(2’^"^l{z[ 


(p) 


> 0 }), and thus. 


G-M 


e = 


27r 


r(i-I)(2 Cr(-y)|co.(f)|) 


where we have also used the expression for E(Z(^) given in Remark 


3.4 


and 


f = _^(M + G)E(4^)4(p) (Z'f^)) + 


7{p)^ 


Gr(-y) 


{M^ -G^). 


(5.4) 


(5.5) 


In the presence of a continuous component. Theorem |4.1| and Corollary |4.3| supply short-term expansions for ATM 
digital call prices and the implied volatility skew. The coefficients dfe, e, and /, are as in (4.3)-(4.5), but in the 
zero-correlation case, the last of these reduces to 




/ = 


-(-!)( 


-MG{l) + GG{-l) 
y -1 


G(l) + G(-l) 1 
a^Y 2 


(5.6) 


We also note that the first two d^-coefRcients are given in Remark 4.2 (f), whereas in the CGMY-case, G(l) = G(—1), 
and all the dkS vanish. 


To assess the accuracy of the approximations, we compare them to the “true” values of the quantities, which are 
estimated using Monte Carlo simulation. First, for ATM digital call options, we use the measure transformation 
introduced in Section [2] to write 


P(Yt + Ft > 0) = E(e-^‘l{x,+v*>o}) = E(e 


-Mz4+Gz)"7 






(5.7) 


where Z^f^ and are strictly Y-stable random variables with Levy measures of the form (2.10 ^ Such variables can 
be simulated efficiently, and standard discretization schemes can be used for the continuous component F, so (5.7) can 
be used to obtain an unbiased estimate of ATM digital call prices. For the implied volatility slope, we first estimate 
near-the-money option prices and implied volatility, using a Monte Carlo procedure similar to the one based on (5.7) 
for digital options, and then use those values to numerically estimate the ATM volatility slope. 


As mentioned in the introduction, a common drawback of short-term expressions is that their convergence may be 
slow, and only satisfactory at extremely small time scales. However, the performance is highly parameter-dependent, 
and in this section we consider parameter values that are of relevance in financial applications, based on results in [SJ 
Tables 1 and 5], and [33 p.82], where the tempered stable model is calibrated to observed option prices. 

Figure [^compares the asymptotic expansions for ATM digital call prices, under a pure-jump tempered stable model, 
to the true price estimated by Monte Carlo simulation. The axes are on a logj^Q-scale, with time-to-maturity in years. 
The first and second order approximations are of order and f, respectively, and it is clear that the second order 

approximation significantly improves the first order approximation, and gives good estimates for maturities up to at 
least one month. Figure then shows the implied volatility smile and the ATM slope approximation. The maturity is 
0.1 years (over one month), and the slope approximation captures the sign of the slope, and even gives a good estimate 
of its magnitude: In Panel (a), the slope approximation is 0.320 compared to a Monte Carlo estimate of 0.328; in Panel 
(b) the approximation gives —0.456 while the Monte Carlo estimate is —0.521. 


^Equivalently, and are F-stable random variables with location parameter 0, skewness parameters 1 and —1, and scale parameters 
(tC(l)|cos(7ry/2)|r(-y))^''^ and (tC(-l)| cos(7ry/2)|r(-y))^/’^, respectively. 
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log(t) 


Figure 1. ATM digital call option prices computed by Monte Carlo, and the first- and second-order ap¬ 
proximations, under a pure-jump tempered stable model. Time is in years and both axes on a logjQ-scale. 
(a) ((7(1), C(-l), G, M, y) = (0.0088,0.0044,0.41,1.93,1.5) as suggested in [5]. (b) (G(l), G(-l), G, M, T) = 
(0.015,0.041,2.318,4.025,1.35) as suggested in [25]. 




Figure 2. The volatility smile (red) as a function of log-moneyness, and the second-order slope approximation 
(blue). Time-to-maturity is 0.1 years, and the models are the same as in Figure]^ 


Figure [^carries out the same analysis for a tempered stable model with a nonzero Brownian component. Panel (a) 
compares the Monte Carlo estimates of ATM digital call prices to the first and second order approximations, which 
are of order and respectively. As in the pure-jump case, the second order approximation gives a good 

estimate for maturities up to at least a month. Panel (b) then displays (in red) the volatility smile and the ATM 
slope approximation for maturity t = 0.1 years (over 1 month). For comparison, it also shows (in blue) the volatility 
smile and slope approximation when the constant volatility Brownian component is replaced by a Heston stochastic 
volatility process with the same spot volatility, and leverage parameter p = —0.3. In both cases the approximations 
clearly capture the sign of the slope, and are in fact very close to the Monte Carlo estimates: 0.112 (resp. 0.305) 
compared to Monte Carlo estimates of 0.116 (resp. 0.289) in the Brownian case (resp. Heston case). 

Finally, for completeness, we include in Figure|^a comparison of the first and second order approximations for ATM 
call option prices from m, and the true values estimated by Monte Carlo simulation. 

5.2. Empirical application: S&:P500 implied volatility skew. In this section we analyze the implied volatility 
skew in S&P500 options, and compare it to the results of Sections and Our dataset consists of daily closing bid 
and ask prices for S&P500 index options, across all strikes, K, and maturities, t, from January 2, 2014, to January 
31, 2014 (21 business days). We take the mid-values of bid and ask prices as our raw data, and observations with 
time-to-maturity of less than five days are dropped to minimize the impact of microstructure effects. For each day and 
maturity we also visualize the quoted prices to check for obvious outliers. 

We base the construction of the implied volatility curve on out-of-the-money (OTM) options since in-the-money 
(ITM) options are infrequently traded compared to OTM options and, thus, their prices are typically less reliable. 
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(a) 


X 

Q. -0.297 


- Monte Carlo Price 

- 1st-order Approx. 

- ^nd-order Approx. 


log(t) 



Figure 3. (a) ATM digital call option prices computed by Monte Carlo, and the first- and second-order 
approximations. Time is in years and both axes on a logj^g-scale. The model is tempered stable with a Brownian 
component and parameters ((7(1), C(—1), G, M, V, ct) = (0.0040,0.0013,0.41,1.93,1.5,0.1) as suggested in [5]. 
(b) The red curves show the volatility smile and the second-order slope approximation for maturity 0.1 years. 
The blue curves show the same quantities after replacing the Brownian component by a Heston process with 
the same spot volatility and p = — 0.3. 
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Figure 4. ATM call option prices computed by Monte Carlo, and the first- and second-order approximations 
from m- Panels (a) and (b) correspond to the pure-jump cases in Figures [l|^ while panel (c) corresponds to 
the mixed case in Figure]^ 


More specifically, we follow a standard procedure, described in [3] and |42j . in which the put-call parity and liquid 
at-the-money (ATM) options are first used to compute the implied forward price of the underlying Then, OTM 
options are used to compute the implied volatility for different strike prices. That is, put options (resp. call options) 
are used for strike prices that are below (resp. above) the forward price. 

Figure shows some stylized features of our data. Panel (a) shows how the strike prices of the typically liquid 
25-delta option^ become increasingly concentrated around the ATM strike, as time-to-maturity decreases, which is 
one reason for the importance of considering a small-moneyness regime in short-time. Panel (b) then shows how the 
implied volatility smile becomes increasingly skewed as time-to-maturity decreases. It also clearly shows that the left 
wing (corresponding to OTM put options) is steeper than the right wing (corresponding to OTM call options), which 
has consistently been observed in S&P option prices since the market crash of 1987, and reflects the negative skewness 
in the underlying distribution of risk-neutral returns, or, equivalently, the high demand for protective put options 
against downward index movements. 


^The ATM strike is taken to be the strike price at which the call and put options prices are closest in value. We also set the risk-free 
interest rate to zero, but using a nonzero rate based on U.S. treasury yields did not change the results of our analysis since the rate is close 
to zero over the sample period and the time-to-maturity is small. 

®The 25-delta put (resp., call) is the put (resp., call) whose strike price has been chosen such that the option’s delta is -25% (resp., 25%). 
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As explained in previous sections, the short-term behavior of the skew differs significantly from one model setting to 
the next: In purely continuous models the skew is bounded, but in jump-models the relationship between the skew and 
time-to-maturity is a power-law, where the exponent depends on properties of the jump-component. This observation 
invites a model selection and calibration procedure based on a comparing the skew observed in real markets, to the 
model-skew under different model assumptions. This is similar in spirit to the approach in El, where the short-term 
decay of option prices is used to infer properties of the underlying asset price process. 


First, we need to decide on a measure of the implied volatility skew. Various skew-measures have been proposed 
in the literature (see, e.g., [30]). In this work, we simply estimate the ATM skew by taking the slope of the volatility 
smile between two near-the-money options. For consistency we alwa ys c hoose OTM 25-delta options, since they are 
actively tradec^ and are indeed near-the-money options (see Figure I.l|. Also, in order to fairly compare the skew 
across different dates and overall volatility levels, all volatilities for a given date are normalized by the CBOE Volatility 
Index (VIX), which is a measure of the average volatility of the S&P index. 


Figure l^a) shows the estimated skew for options with maturities up to one year, for each day in our dataset. The 
power-law behavior is evident, not only for short-term options, but it seems to continue to hold for longer maturities 
as well. Consequently, in Figure [^b) the log-skew seems to be linear in log-maturity. Thus, to estimate the exponent 
of the power law, we run a linear regression of log-skew on log-maturity, and take the slope coefficient as an estimate of 
the exponent of the power law. Our results can then be used for model selection and calibration: A nonnegative slope 
coefficient implies a purely continuous model or a jump-diffusion model, while a negative slope coefficient is consistent 
with a jump component of infinite activity. In the latter case, the magnitude of the slope coefficient can also be used 
to distinguish between jump-components of finite and infinite variation. 


We carry out this regression for each of the 21 days in our dataset, using options with time-to-maturity less than 
0.25 years (3 months), and days with at least four maturities less than 0.25 years (which discards only one day). We 
find that the linear model fits the data extremely well, with the average i?-squared being 0.98. Moreover, the slope 
coefficient is consistently observed to be between —0.3 and —0.4, with an average of -0.36, which, as explained above, 
contradicts a jump component with finite activity, as well as an infinite activity jump-component of finite variation, in 
which case the order of the skew is -0.5 (cf. [15]). On the other hand, this slope coefficient is in line with the results for 
the infinite variation jump-models studied in Sections m and can be used to calibrate important model parameters. 


In particular, this yields a simple procedure to calibrate the index of jump activity of the process, V, which can 
be viewed as a new forward-looking tool to assess this fundamental parameter, complementing the popular rear-facing 
estimation methods based on high-frequency observations of the underlying asset’s returns (cf. [1]). First, in the pure- 
jump case of Section]^ the order of the skew is —1/2 if (7(1) ^ C{—1), but 1/2 — 1/V S (—1/2,0) otherwise, so our 
regression results suggests that (7(1) = (7(—1) and Y £ (1.11,1.25). However, for the mixed model of Section]^ the 
skew is bounded if (7(1) = (7(—1), but of order 1/2 — Y/2 G (—1/2,0) otherwise, so the regression results, together 
with the sign of the skew (see Remark 4.4), suggest that (7(1) < (7(—1) and Y G (1.6,1.8). These results show that 
both the pure-jump model and the mixed model are able to capture the short-term order of the skew, but the latter 
value of Y is perhaps of greater interest, as several studies point to the presence of a continuous component in the 
returns process, and Y > 1.5 for actively traded stocks (cf. U, [2], [TT]j. 


In the analysis above we used options with time-to-maturity less than 3 months. However, one of the main conclusion 
of the first part of this section was that the asymptotic expansions seem to give good approximations for options with 
maturities up to 0.1 years (just over one month). We therefore repeat the analysis using only options with maturities 
below or around the one month mark, excluding days with fewer than three maturities in that range, which eliminates 
3 days out of 21. In this case, the average i?-squared is 0.97, and the average slope is -0.31, which, in the pure-jump 
case, indicates (7(1) = (7(—1) and Y = 1.23, while, in the mixed-case, it again suggests (7(1) < (7(—1) and Y = 1.62. 
Those Y -values are of similar magnitude as in the previous analysis, but slightly more moderate in the sense that they 
are closer to 1.5. We also note that the analysis based on options with maturities shorter than 0.1 years resulted in 
more stable regression results between days, compared to when including maturities shorter than 0.25, i.e. most of the 
slope coefficients were quite close to the average -0.31. 


Finally, let us remark that even though one cannot reject one of our models in favor of the other based solely on 
the behavior of the skew, they can be distinguished, in principle, by the fact that the ATM volatility converges to 
zero in the pure-jump case, but to the nonzero spot volatility in the mixed case. As already described in El, this 
property of models with a continuous component can also be exploited when using short-term expansions to calibrate 


^Repeating the analysis using 10-delta options did not have a qualitative effect on the outcome. 
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Figure 5. (a) The log-moneyness (k = \n{K/F), where F is the forward price) of outstanding 25-delta call 
options (k > 0) and put options (n < 0) in Jan 2014. (b) The implied volatility smiles on Jan 15, 2014, 
corresponding to maturities ranging from 0.012 (3 days) to 0.25 (3 months). 




Figure 6. (a) The ATM implied volatility skew as a function of time-to-maturity, for each business day in 
Jan 2014. (b) The natural logarithm of the skew in panel (a) as a function of the natural logarithm of 

t ime-t o- mat urity. 


their parameters, to circumvent the fact that the spot volatility is not directly observable (see, e.g., [29]). Figure 
l^a) displays for each day the ATM volatility as a function of time-to-maturity, and a quick inspection indicates that 
extrapolating to the zero-maturity does not seem to result in zero volatility, as would be the case in a pure-jump 
model. Furthermore, Figure [^b) shows how the ATM implied volatility for the shortest outstanding maturity moves 
in tandem with the corresponding VIX measurement. 


6. Conclusions and future work 

As time-to-maturity becomes short, it is observed empirically that the liquid strike prices become increasingly 
concentrated around the ATM strike. The short-term volatility smile is therefore sometimes analyzed in terms of three 
quantities: the ATM implied volatility level, skew, and convexity (see, e.g., m and |12|). 

The present work focuses on the ATM implied volatility skew (i.e. the strike-derivative). We obtain high-order short¬ 
term expansions for the skew under models with a Levy jump-component of infinite variation, and in the presence of a 
nonhomogeneous continuous component we quantify explicitly the short-term skew-effects of both jumps and stochastic 
volatility. Our proofs utilize a relationship between the skew and transition probabilities of the form P(S't > Sq), i.e. 
prices of ATM digital call options, and as auxiliary results we also obtain short-term approximations for the delta 
of ATM options. Simulation results indicate the validity of our results for options with maturities up to at least one 
month, and we show that the volatility skew in recent S&P500 index options is in accordance with the infinite variation 
jump-component of our models. 
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Figure 7. (a) ATM implied volatility as a function of time-to-maturity for each business day in Jan 2014. 
(b) ATM implied volatility for the shortest outstanding maturity compared to the VIX index. 


It is also natural to wonder about the short-time behavior of the ATM smile convexity, which, just like the skew, 
is of great importance in financial markets (see Section O- Specifically, the ATM convexity is defined as the second 
order strike-derivative of the implied volatility. 


d'^cT {k, t) 




( 6 . 1 ) 




and, following similar steps as those used to derive (1.3), it is possible to show that the key quantity needed to analyze 


the convexity is /sj(S'o), where fs^ is the risk-neutral probability density of St- Moreover, just like the transition 
probability F{St > Sq) could be linked to the delta of ATM options, the probability density fsti^o) is connected 
to the gamma of ATM options. In order to analyze the behavior of fsti^o), it seems necessary to also analyze the 
characteristic function of the log-returns process and, more specifically, its decay properties. This is in sharp contrast 
to the methods used in the present work that do not rely on inverse Fourier representations of the option prices and 
transition probabilities. The latter approach is left for future work. 


Appendix A. Additional Proofs 


Lemma A.l. Let V be as in (1.10)), with n {Yt) and cr{Yt) replaced by and at, defined in {f.lf )■ Also let a't, a'(, 
at, and y*, be the stopped processes in (4-30), and af := \J\4lalds. Then the following relations hold for any p > 1: 


(i) E \flt - fiof = 0{ti), t^O. 
{ii) E \at - aof = 0{t^), t ^ 0. 
(hi) E\af — aof = 0{t^), t —?> 0 . 



(v) For^l,^f, and , as in (4.30)-(4.31), wehaveE\f4\ = 0(t) and E |^^|-|-E|^( — = 0{t^), t —>■ 0. 


Proof. Let L be a common Lipschitz constant for p,t, at, and 74 . 

(i) By the Lipschitz continuity of fit, and the Burkholder-Davis-Gundy (BDG) inequality, we can find a constant Cp 
such that 

E I/I 4 - Mor < LPE - 2 /or < LPCp{E(^j\sdsy +E(^j\^^ds) = 0 (tS), t ^ 0, 

since dg and 7 s are bounded. 

(ii) is proved in a similar way, and for (hi) we use the boundedness of at, Jensen’s inequality, and (ii) to write 

® ■ "“I" - I *)' ^ A)'] I ^ 0- 
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(iv) We can write E((Tf — ao) = E((ctJ‘)^ ~ / (2o’o) +®^(((^t ) ~ + cro) — l/(2o'o))), where the second term 

is 0{t) by (iii), while for the first term we have by Ito’s lemma 

rt 


E(( 


CTf - CTn = 




+ 


'^dujds'j = 0{t), t —>■ 0, 


due to the fact that the expected value of the stochastic integral is zero, 
(v) By Cauchy’s inequality and Ito’s isometry we have 


E < 


Similarly, 




\lI e(^ {(ruOiu + ^cr'uli)duJ ds = 0{t^, t 

J a'^joYduds = 0{ti), t- 


0 . 


0 , 


because by the boundedness of cr(j and 'ju, we can find a constant K such that 

E {d-uJu - o-o7o)^ < KE - "tof + KE (ct^ - < 2LKE (Wat - Vof = 0{u), u-)■ 0, 

where in the last step we again used the BDG inequality. Similarly, Cauchy’s inequality and Ito’s isometry yield 

E|Ct| =0(0- 

(vi) follows from the triangle inequality and the following three identities. First, by (iii) above, we have 


E 


- O-Q - 


K) 


2(To 


= ^E(at - dof = 0(t), t 

2(70 


0 . 


Second, by Ito’s Lemma, 


E 


K)^-c^o-7/ [ ‘^(^uCr'uludWlds <\ [ [ E 2cr„cr>„ + ^ ((cr(,)^ + ct^ct") 7^ duds = 0(t), t ^ 0, 
t Jo Jo ^ Jo Jo 2 


since the integrand in the last integral is bounded. Third, Cauchy’s inequality and Ito’s isometry can be used to show 


E 


/ / {o'ufTu'lu - o-oo-'olo) dW^ds = 0{t), t 0, 

Jo Jo 


by following similar steps as in the proof of (v). 


□ 


References 

[1] Y. Ait-Sahalia and J. Jacod. Estimating the degree of activity of jumps in high-frequency data. Annals of Statistics, 37(5A), 2202—2244, 
2009. 

[2] Y. Ai't-Sahalia and J. Jacod. Is Brownian motion necessary to model high-frequency data? Annals of Statistics, 38:3093-3128, 2010. 

[3] Y. Ai't-Sahalia and A. Lo. Nonparametric estimation of state-price densities implicit in financial asset prices. Journal of Finance, 53(2), 
1998. 

[4] E. Alos, J. Leon, and J. "Vives. On the short-time behavior of the implied volatility for jump-diffusion models with stochastic volatility. 
Finance and Stochastics, 11(4), 571-589, 2007. 

[5] L. Andersen and A. Lipton. Asymptotics for exponential Levy processes and their volatility smile: Survey and new results. Int. J. Theor. 
Appl. Finance, 16(1), 1-98, 2013. 

[6] D.S. Bates. The crash of ’87: Was it expected? The evidence from options markets. Journal of Finance, 46(3), 1009-1044, 1991. 

[7] G. Bakshi, N. Kapadia, and D. Madan. Stock return characteristics, skew laws, and the differential pricing of individual equity options. 
Review of Financial Studies, 16(1), 101-143, 2003, 

[8] L. Bergomi. Smile Dynamics. Risk magazine, 9, 117-123, 2004. 

[9] J. Bertoin. Levy processes. Cambridge University Press, 1998. 

[10] P. Carr, H. Geman, D. Madan and M. Yor. The fine structure of asset returns: An empirical investigation. Journal of Business, 75(2), 
303-325, 2002. 

[11] P. Carr and L. Wu. What type of process underlies options? A simple robust test. Journal of Finance, LVIII(6), 2003. 

[12] R. Cont and P. Tankov. Financial modelling with jump processes. Chapman Sz Hall, 2004. 

[13] L. De Leo, V. Vargas, S. Ciliberti and J.-P. Bouchaud. We’ve walked a million miles for one of these smiles. Preprint, available at 
http://arxiv.org/abs/1203.5703, 2012. 

[14] G. Fusai and A. Meucci. Pricing discretely monitored Asian options under Levy processes. Journal of Banking & Finance, 32, 2076- 
2088, 2008. 

[15] S. Gerhold, LC. Giilum, and A. Pinter. The small-maturity implied volatility slope for Levy models. Available at arXiv:1310.3061 
2014. 













SKEW ASYMPTOTICS UNDER STOCHASTIC VOLATILITY AND LEVY JUMPS 


29 


[16] J. Fajardo and E. Mordecki. Skewness premium with Levy processes, Quantitative Finance, 14(9), 1619—1626, 2014. 

[17] J.E. Figueroa-Lopez and M. Forde. The small-maturity smile for exponential Levy models, SIAM Journal on Financial Mathematics, 
3(1), 33-65, 2012. 

[18] J.E. Figueroa-Lopez, R. Gong, and C. Houdre. High-order short-time expansions for ATM option prices of exponential Levy models. 
To appear in Mathematical Finance, DOI: 10.1111/mafi.12064, 2014. 

[19] J.E. Figueroa-Lopez, R. Gong, and G. Houdre. Small-time expansions of the distributions, densities, and option prices under stochastic 
volatility models with Levy jumps. Stochastic Processes and their Applications, 122, 1808-1839, 2012. 

[20] J.E. Figueroa-Lopez and G. Houdre. Small-time expansions for the transition distribution of Levy processes. Stochastic Processes and 
their applications, 119, 3862—3889, 2009. 

[21] J.E. Figueroa-Lopez and S. Olafsson. Short-time expansions for close-to-the-money options under a Levy jump model with stochastic 
volatility. To appear in Finance and Stochastics, DOI: 10.1007/s00780-015-0281-z, 2015. 

[22] K. Gao and R. Lee. Asymptotics of implied volatility to arbitrary order, Finance and Stochastics, 18(2), 349—392, 2014. 

[23] J. Gatheral. The volatility surface: A practitioner's guide. Wiley finance series, 2006. 

[24] O. Kallenberg. Foundations of modern probability. Springer-Verlag, Berlin, New York, Heidelberg, 1997. 

[25] R. Kawai. On sequential calibration for an asset price model with piecewise Levy processes. lAENG International Journal of Applied 
Mathematics, 40(4), 239-246, 2010. 

[26] M. Konikov and D. Madan. Stochastic volatility via Markov chains. Review of Derivatives Research, 5, 81-115, 2002. 

[27] 1. Koponen. Analytic approach to the problem of convergence of truncated Levy flights towards the Gaussian stochastic process. 
Physical Review E, 52, 1197-1199, 1995. 

[28] R. Lee. Implied volatility: Statics, dynamics, and probabilistic interpretation. Recent advances in applied probability, 241—268, Springer, 
New York, 2005. 

[29] A. Medvedev and O. Scailllet. Approximation and calibration of short-term implied volatility under jump-diffusion stochastic volatility. 
Review of Financial Studies, 20(2), 427-459, 2007. 

[30] S. Mixon. What Does Implied Volatility Skew Measure? Journal of Derivatives, 18(4), 9—25, 2011. 

[31] A. Mijatovic and P. Tankov. A new look at short-term implied volatility in asset price models with jumps. To appear in Mathematical 
Finance, Available at arXiv: 1207.0843, 2013. 

[32] J. Muhle-Karbe and M. Nutz. Small-time asymptotics of option prices and first absolute moments. Journal of Applied Probability, 
48(4), 1003-1020, 2011. 

[33] S. Olafsson, Applications of short-time asymptotic methods to option pricing and change-point detection for Levy processes. Ph.D. 
thesis, Purdue University, 2015. 

[34] J. Pan. The jump-risk premia implicit in options: Evidence from an integrated time-series study. Journal of Financial Economics, 63, 
3-50, 2002. 

[35] M. Roper and M. Rutkowski. On the relationship between the call price surface and the implied volatility surface close to expiry. Int. 
J. Theor. Appl. Finance, 12, 427-441, 2009. 

[36] M. Rosenbaum and P. Tankov. Asymptotic results for time-changed Levy processes sampled at hitting times. Stochastic processes and 
their applications, 121, 1607—1633, 2011. 

[37] K. Sato. Levy processes and infinitely divisible distributions. Cambridge University Press, 1999. 

[38] W. Schoutens. Levy processes in finance. Wiley, 2003. 

[39] P. Tankov. Pricing and hedging in exponential Levy models: Review of recent results. Paris-Princeton Lecture Notes in Mathematical 
Finance, Springer, 2010. 

[40] Y. Xing, X. Zhang, and R. Zhao. What does individual option volatility smirk tell us about future equity returns? Journal of Financial 
and Quantitative Analysis, 45(3), 641—662, 2010. 

[41] S. Yan. Jump risk, stock returns, and slope of implied volatility smile. Journal of Financial Economics, 99(1), 216-233, 2011. 

[42] J.E. Zhang and Y. Xiang. The implied volatility smirk. Quantitative Finance, 8(3), 263-284, 2008. 

[43] V. M. Zolotarev. One-dimensional stable distributions. Amer. Math. Soc., Providence, R.L, 1996. 


Department of Mathematics, Washington University in St. Louis, St. Louis, MO 63130, USA 
E-mail address: figueroa@inath.wustl.edu 

Department of Statistics and Applied Probability, University of California, Santa Barbara, CA 93106, USA 
E-mail address: olafsson@pstat.ucsb.edu 


